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CHAPTER  6 


TASK  4:  PROTOTYPE  DESIGN;  REAL  TIME  BREADBOARD  SIMULATION 


This  volume  Is  the  second  of  four  which  collectively  summarize 
the  results  from  the  development  of  a  control  strategy  for  ejection 
seats  for  the  Air  Force  Program  entitled  "Vectored  Thrust  Digital 
Flight  Control  for  Crew  Escape."  Volume  I  of  this  effort  concentrated 
on  posing  the  fundamental  control  problem  and  reviewed  several  control 
synthesis  techniques  for  potential  application  to  this  problem.  The 
"acceleration  control"  approach  based  on  Meyer  (1975)  best  met  the 
stringent  time  requirements  In  dealing  directly  with  the 
life-threatening  accelerations  typical  in  the  high  dynamic  environment 
of  ejection  seat  deployment.  Subsequent  discussion  dealt  with 
adapting  Meyer's  approach  to  the  specific  requirements  of  the  ejection 
seat  problem  and  detailed  the  benefits  of  the  selected  control  design. 
Also  discussed  in  Volume  I  was  the  result  of  a  technology  survey  that 
Investigated  the  current  availability  of  sensors,  actuators,  and 
microprocessors  to  meet  the  respective  state  estimation  accuracy, 
control  energy  and  computational  requirements  Inherent  in  the 
implementation  of  the  new  control  design. 

Volume  II  focuses  on  the  description  of  the  simulation  effort  and 
results  obtained  using  the  "acceleration  control"  approach  applied  to 
the  ACES  II  ejection  seat  given  the  confines  of  existing  technologies. 
The  primary  goal  of  Task  4  is  the  evaluation  of  the  feasibility  of 
real-time  operation  of  a  microprocessor  based  implementation  of  the 
control  design  when  applied  to  a  wide  variety  of  escape  conditions 
ranging  from  benign  to  immediately  life  threatening. 


Given  the  enoraous  computaclonai  overhead  in  generating  the  truth 
model  for  high  fidelity  flight  simulations,  standard  digital  computers 
are  typically  unsuitable  for  real-time  operation.  That  is,  the  truth 
models  require  generation  of  the  necessary  sensor  Inputs  for  use  by 
the  control  law  and  the  evaluation  of  the  subsequent  airframe  response 
to  control  outputs.  For  high  dimensional  models  operating  at  fast 
rates  (on  the  order  of  SO  to  100  Ht), the  requirements  would  exceed 
the  capacity  of  most  digital  computers  with  the  exception  of  large 
dedicated  main  frames. 

The  approach  taken  here  is  to  simulate  in  analog  hardware  as  many 
of  the  models  as  practicable  with  the  exception  of  the  digital  control 
law  under  test.  The  analog  computer  is  by  nature  a  parallel 
processing  device  so  that  real-time  simulation  of  even  enormously 
complex  system  models  is  feasible.  In  this  "hybrid”  mode,  the 
microprocessor  based  system  solicits  sensor  inputs  directly  over  an 
analog  to  digital  converter  (M)C)  interface  and  delivers  the  quantized 
control  outputs  to  the  actuation  system  via  digital  to  analog 
converters  (OAC).  The  interfaces  then  accurately  reflect  the  actual 
communication  medium  of  the  operational  host  system.  That  is,  the 
effect  of  signal  quantization,  sampling  error,  ambient  electrical 
noise,  hardware  response  time  and  control  cycle  timing  are  all 
implicit  in  system  operation. 

The  simulation  test  cases  examined  are  meant  to  exercise  the 
control  law  over  the  entire  spectrum  of  flight  conditions  for  which  it 
is  designed.  In  the  low  altitude  regime  the  flight  scenarios  exercise 
"adverse  attitude"  conditions,  l.e.,  conditions  In  which  ground 
collision  is  likely  If  the  opportunity  for  immediate  corrective  action 
is  delayed.  High  dynamic  pressure  operation  is  demonstrated 
in  low  altitude  and  high  altitude  situations  while  other  test 
cases  are  aimed  at  examining  the  sensitivity  of  control  performance  to 
large  angle  variations  In  Initial  attitude.  Of  paramount  Importance 
in  the  evaluations  Is  the  control  system  effectiveness  In  maintaining 


Che  "acceleraClon  radical"  below  the  lethal  or  high  probability  of 
Injury  Halt  while  demonstrating  proper  terminal  attitude  control  for 
subsequent  phases  of  chute  deployment  and  recovery. 

Volume  II  is  subdivided  into  three  chapters,  6  through  8.  Chapter 
6  is  devoted  to  the  description  of  the  hybrid  simulation  models  and 
Includes  a  discussion  of  Che  aerodynamics,  dynamics  equations,  plus 
sensor  and  actuator  models  comprising  the  analog  segment  of  the 
simulation  (6.2).  The  digital  segment  representing  Che  control  and 
actuator  element  logic  is  discussed  in  section  6.3.  The  allocation  of 
hardware  to  perform  the  hybrid  processing  is  described  in  section  6.4 
and  includes  a  brief  description  of  Che  Martin  Marietta  computing 
configuration  as  well  as  Che  microprocessor  based  system  developed  by 
Unidynamics.  Overlap  with  some  material  presented  in  other  volumes  of 
this  report  Is  inevitable  but  for  Che  sake  of  continuity  and 
completeness  some  information  is  repeated  when  necessary. 

Chapter  7  represents  the  main  body  of  the  results.  The  purpose 
of  Chapter  7  is  to  demonstrate  that  the  major  control  objectives  are 
met  when  evaluated  in  "difficult"  situations.  Aside  from  the  emphasis 
on  control  system  performance,  ocher  results  are  presented  Chat 
demonstrate  consistency  between  alternate  analog 
simulation  models.  For  example,  the  hybrid  simulation  allows  for 
generation  of  results  through  the  use  of  a  main  frame  computer 
replacing  the  microprocessor  Implementation  of  the  control  algorithm. 
Close  agreement  in  the  results  is  obtained  with  the  alternate 
configurations  as  demonstrated  in  section  7.2.  Also  Included  as  a 
comparison  in  the  results  are  the  effects  of  "open  loop"  operation  of 
Che  system.  Here  "open  loop"  refers  to  fixed  force  magnitude  of  the 
main  thruster  with  no  means  of  attitude  control  corresponding  to 
current  ACES  II  performance  due  to  the  omission  of  the  STAPAC  and  the 
drogue  from  the  system.  These  comparisons  vividly  demonstrate  the 
direct  benefits  of  intelligent  ejection  seat  control. 


The  results  exercising  the  MZL'SPEC  cases  as  defined  In  are 
presented  la  section  7.3  .  Those  cases  reflect  the  "adverse  attitude" 
cases  previously  aentloned.  Section  7.4  is  devoted  to  examining  the 
alternate  20  cases  as  required  by  SOW  paragraph  4. 5. 4. 2  which 
exercises  the  control  law  over  a  wide  range  of  Initial  speed,  altitude 
and  attitude  conditions. 

Paragraph  7.6  Is  a  discussion  of  the  memory  and  throughput 
measurements  for  the  Motorola  68000  microprocessor  (MC  68000)  In  this 
control  application.  The  original  Intent  of  Implementing  the  control 
algorithm  In  a  higher  order  language  (PASCAL)  has  proven  to  yield 
results  far  from  satisfactory  with  respect  to  system  timing.  However 
an  assembly  version  of  the  control  code  has  proven  successful  In 
meeting  the  stringent  timing  requirements  of  the  duty  cycle. 

Paragraph  7.7  analyzes  the  stability  of  the  control  system  to 
variations  In  parameters  such  as  gains,  sensor  and  actuator 
bandwldths,  as  well  as  sensitivity  to  Initial  conditions  and 
uncertainty  in  moment  of  Inertia  and  CG  location. 

Finally,  paragraph  7.8  la  a  stand  alone  section  which  evaluates 
the  effects  of  pilot  CG  motion  on  control  system  performance.  Due  to 
limitations  In  the  pilot  harness  restraints,  some  limited  pilot  motion 
when  subjected  to  high  g's  Is  Inevitable.  As  a  consequence  the  net 
seat/pilot  center  of  gravity  and  moment  of  Inertia  is  actually  a  time 
varying  quantity.  The  control  law  maintains  a  constant  estimate  of 
these  key  system  parameters  so  that  the  variations  In  performance  Is 
of  some  Issue  In  determining  robustness  of  the  control  approach.  The 
results  Included  In  paragraph  7.8  demonstate  an  Insignificant  level  of 
sensitivity  to  pilot  motion  disturbance. 

Chapter  8  offers  some  basic  conclusions  on  the  overall  system 
performance  plus  Indicated  restrictions  for  deployment  and  offers 
recommendations  for  further  studies  of  and  enhancements  in  the  control 


system  design 


6, 1  Overview  of  the  hybrid  simulation  model 

The  hybrid  simulation  for  the  ejection  seat  problem  Is  a 
collection  of  modules  that  model  the  dynamics  and  response  for  the 
components  ot  the  ACES  II  ejection  seat  flight  system.  Three 
classifications  of  models  are  evident:  1)  environment  models,  2) 
sensor  and  actuation  system  models  and  3)  the  digital  control 
algorithm.  The  environment  models  are  concerned  with  representing 
with  the  highest  degree  of  realism  the  effects  of  the  aerodynamic 
forces  and  torques  acting  on  the  ejection  seat  as  a  function  of  the 
specific  flight  conditions,  l.e.  ,  ejection  seat  altitude,  speed  and 
attitude.  The  sensor  models  receive  the  best  estimate  of  the  "truth" 
from  the  ejection  seat  dynamics  model  and  in  turn  generate  realistic 
signals  which  incorporate  the  major  sources  of  instrument  error.  Here 
the  most  pressing  source  of  error  is  instrument  delay  since  the 
control  system  response  time  is  critical  in  alleviating  the 
life-threatening  forces  acting  on  the  pilot's  bodv.  The  actuation 
system  models  are  similarly  concerned  with  imposing  the  natural 
instrument  delay  in  responding  to  control  law  command  inputs.  Hence, 
the  actual  forces  and  torques  sensed  by  the  seat  differ  from  the 
idealized  commands  due  to  the  limitations  of  the  actuation  system. 

The  digital  control  law  accepts  the  raw  Instrument  inputs  and,  given 
certain  control  parameters  (gains,  time  constants,  trajectory  profile 
specifications),  generates  the  idealized  control  signals  for  use  by 
the  actuation  system  just  mentioned. 

In  digital  simulations  of  continuous  time  systems  the  model 
components  are  discrete  in  form,  i.e. ,  at  each  discrete  time  step  the 
individual  modules  evaluate  an  average  value  for  projection  over  the 
ensuing  time  step  which  are  mathematically  combined  to  propagate  the 
system  states.  In  the  hybrid  simulation  the  approach  is  to  reduce  the 
mathematical  model  to  equivalent  analog  circuitry  and  as  a  result 


represent  the  systen  states  continuously  In  tlne«  The  digital  control 
law  under  teat  Is  by  definition  a  discrete  evaluation  entity  so  that 
the  continuous  system  states  are  sampled  at  fixed  Intervals  In  time 
and  the  associated  quantized  control  outputs  generated  for  subsequent 
use  in  state  propagation  by  continuous  time  models. 

The  Integrated  model  components  of  the  hybrid  ejection  seat 
simulation  and  control  system  Is  Illustrated  In  Figure  6.1  while  the 
associated  parameters  of  the  models  as  well  as  chose  in  the  following 
paragraphs  are  defined  In  Table  6. L.  A  qualitative  description  of  the 
system  In  Figure  6.1  is  warranted  to  clarify  Che  necessity  of  Che 
Illustrated  model  components. 

At  Che  heart  of  the  integrated  simulation  Is  the  six  degree  of 
freedom  flight  model  (6D0F)  which  serves  as  Che  "truth"  model  for 
evolution  of  the  key  system  states.  With  the  initial  conditions  for 
Che  ejection  seat  defined,  Che  6D0F  flight  model  integrates  the 
combined  aerodynamic  forces  and  torques  with  Che  achieved  actuator 
forces  and  torques  to  update  the  "truth"  values  of  the  vehicle  dynamic 
states  (u,v,w,x,y,z,q,9,:i/,p,qtr).  The  truth  values  of  the  states  are 
represented  in  body  coordinates  ($  ,e  ,4,  ,a^  ,v^  ,yb  ,,'b  j 
processed  by  sensor  models  ,G^^)  which  incorporate 

instrument  error  sources  to  form  Che  final  sensor  estimates 
.V?  ,  )  for  use  by  the  control  algorithm. 

The  control  algorithm  accepts  the  Instrument  estimates  of  the 
vehicle  dynamic  states  plus  specific  control  law  parameters  and 
evaluates  Che  Ideal  desired  body  forces  and  torques  (f^.'^c)  necessary 

for  seat  control.  These  idealized  commands  (fl.T^)  are  transformed 

c 

into  Che  desired  Individual  actuation  system  element  commands 
-1  B  B 

(u  ■  P  lf^,T°))  which  for  the  particular  syj^tem  under  investigation 
are  specifically  rocket  force  magnitudes  (f/.^)  and  the  associated 

Ni  Ni 

pointing  angles  (  9  ,  41  T  which  drive  the  actual  outputs  of  thrust  jet 
^  ^  -B  B 

control  elements.  The  transformations  from  (r  ,Tc  )  to  control 


TABLE  6.1  Definition  of  Terms  in  Hybrid  Simulation 


Block 

Model 

Symbol 

Definition 

Units 

B1 

Aero 

a 

angle  of  attack 

degrees 

S 

side  slip  angle 

degrees 

Cx 

axial  force  coefficient 

none 

Cy 

side  force  coefficient 

none 

^2 

normal  force  coefficient 

none 

Ct 

rolling  moment  coefficient 

$2 

Cm 

pitching  moment  coefficient 

S2 

yawing  moment  coefficient 

S2 

B2 

6  DOF 

U.UQ 

velocity  projected  onto  the  seat 

X  axis 

fps 

o 

> 

> 

velocity  projectd  onto  the  seat 
y  axis 

fps 

W.WQ 

velocity  projected  onto  the  seat 

fps 

z  axis 

P.PO 

rotation  rate  about  the  seat 

X  axis 

R/S 

q.qo 

rotation  rate  about  Che  seat 
y  axis 

R/S 

r.ro 

rotation  rate  about  cne  seat 
z  axis 

R/S 

X.XQ 

inertial  x  position 

ft 

y.yo 

inertial  y  position 

ft 

X.ZQ 

inertial  z  position 

ft 

>.  Oo 

inertial  roll  angle 

R 

9,  9o 

Inertial  pitch  angle 

R 

'i'o 

inertial  yaw  angle 

R 

Vo 

wind  speed 

fps 

• 

U 

acceleration  for  x  seat  axis  for 

fps2 

rotating  observer 


TABLE  6.1  Definition  of  Terms  in  Hybrid  Simulation  (Continued) 


Block 

Model 

Symbol 

Oe finltlon 

Units 

B2  6 

OOF 

• 

acceleration  for  y  seat  axis  for 
rotating  observer 

fps^ 

• 

w 

acceleration  for  z  seat  axis  for 
rotating  observer 

fps2 

• 

p 

rotational  acceleration  about  the 

X  seat  axis 

R/s2 

• 

rotational  acceleration  about  the 
y  seat  axis 

R/s2 

• 

r 

rotational  acceleration  about  the 
z  seat  axis 

R/s2 

• 

X 

X  component  of  Inertial  velocity 

fps 

• 

y 

y  component  of  Inertial  velocity 

fps 

• 

z 

z  component  of  Inertial  velocity 

fps 

0 

$ 

roll  rate 

R/S 

• 

9 

pitch  rate 

R/S 

• 

V 

yaw  rate 

R/S 

Ix'^y 

moment  of  Inertia  about  the  x,y,z 
seat  axes 

slug-ft 

^xz,^xy»^yz 

cross  product  of  Inertia 

slug-ft 

M 

Mach  number 

none 

q 

dynamic  pressure 

lb/ft2 

s 

seat  reference  cross  section  area 

ft2 

o 


seat  (oass 


slugs 


TABLE  6.1  Definition  of  Terms  in  Hybrid  Simulation  (Continued) 


Block  Model  Symbol  Definition 

B2  6  DOF  '  g  local  gravity  acceleration 

magnitude 

R 

X  component  of  resultant  of  rocket 
force 

R 

fy  y  component  of  resultant  of  rocket 

force 

R 

fj  z  component  of  resultant  of  rocket 

force 

b 

a^  vector  of  inertial  accelerations 

expressed  in  body 

xqq  seat  center  of  gravity  in  seat 

units  (x  component) 

yC(j  seat  center  of  gravity  in  seat 

units  (y  component) 

zqq  seat  center  of  gravity  in  seat 

units  (z  component) 

R 

Tjj  X  component  of  rocket  resultant 

torques 

R 

Ty  y  component  of  rocket  resultant 

torques 

R 

Tj  z  component  of  rocket  resultant 

torques 

CONTROL  estimates  of  pitch,  roll,  yaw 

,b 

a^  estimate  of  Inertial  acceleration 

expressed  in  body 

^b 

Vj  estimate  of  Inertial  velocity 

expressed  in  body 


TABLE  6.1  Definition  of  Terms  in  Hybrid  Simulation  (Continued) 


Block 

B4 


Model  Symbol 

Definition 

Units 

,  b 

CONTROL 

escimate  of  inertial  angular  velo¬ 
city  in  body 

R/S 

.b 

escimate  of  inertial  angular 
acceleration  in  body 

R/s2 

b 

“d 

vector  of  desired  angular  velo¬ 
cities 

R/S 

b 

Vd 

desired  velocity  in  body 

fps 

b 

^d 

desired  acceleration  in  body 

fps^ 

•  3 

vector  of  desired  angular 
accelerations 

R/s2 

•'cpo 

control  gain  for  cross  product 
steering 

1/S2 

control  gain  for  cross  product 
steering 

l/S 

translational  acceleration  control 
gain 

none 

^2 

translational  acceleration  control 
gain 

1/S 

k3 

rotational  acceleration  control 
gain 

none 

k4 

rotational  acceleration  control 
gain 

1/S 

rocket  burn  time  constant 

s 

AC 

simulation  time  step 

s 

table  6.1  Definition  of  Terms  in  Hybrid  Simulation  (Continued) 


Block.  Model  Symbol 
B4  control  I 

tki 

?F 

ap 

B5  Config-  P"^ 

uratlon 

_Nl 

J^2 

-^3 

fc 

'^c 

Nl 

N2 

*'c 

N2 

9c 

Nl 

Uc  -  l?c 


Definition  Units 

moment  of  Inertia  slug-ft^ 

control  time  constant  for  l/S 

acceleration 

final  Inertial  pitch  angle  R 

final  inertial  roll  angle  R 

final  inertial  acceleration  fps^ 

Idealized  control  forces  in  body  lb 

idealized  control  torques  in  body  Ib-ft 

inverse  of  control  gradient  dependent 

matrix 

commanded  force  magnitude  for  lb 

rocket  1 

commanded  force  magnitude  for  lb 

rocket  2 

commanded  force  magnitude  for  lb 

rocket  3 

commanded  pointing  angle  ■]>  for  R 

rocket  1 

commanded  pointing  angle  3  for  R 

rocket  I 

commanded  pointing  angle  >  for  R 

rocket  2 

commanded  pointing  angle  3  for  R 

rocket  2 

N2 

,0^.)  control  vector 
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TABLE  6.1  Definition  of  Terns  of  Hybrid  Simulation  (Continued) 


Block  Model 

B6  ACTUATOR 


B7  RESULTANT 


Symbol 

Definition 

/I 

achieved  pointing  angle  for 
rocket  1 

achieved  li;  pointing  angle  for 
rocket  2 

achieved  9  pointing  angle  for 
rocket  I 

6  ^"2 

achieved  9  pointing  angle  for 
rocket  2 

_Nl 

achieved  rocket  force  magnitude 

f 

for  rocket  1 

J*2 

achieved  rocket  force  magnitude 

f 

for  rocket  2 

_N3 

achieved  rocket  force  magnitude 

f 

for  rocket  3 

Ni  Ml  Ni 

fir^  ,5ry  ,5rj 

displacement  of  rocket  1  from  c.g 
in  seat  units 

N3  N3  N3 

6  Tjj  ,  5  Ty  ,  5  Tj 

displacement  of  rocket  3  from  c.g 
in  seat  units 

N2  N2  N2 

displacement  of  rocket  2  from  c.g 
in  seat  units 

b  b 

•  *Ry  »  ‘Rj 

components  of  resultant  rocket 
forces  in  body 

b  b  b 

^Rx  >^Ry  '^Rj 

components  of  resultant  rocket 
torques  in  body 

b 

Cni 


direction  cosine  matrix  from 
rocket  Ni  to  body 


EJECTION  SEAT  SIHULATION/CONTROI.  SYSTEM 


vector  u  Is  obviously  dependent  upon  the  constraints  Inherent  In  the 
c 

particular  actuation  configuration  employed.  Results  employing 
alternate  actuator  configuration  schemes  are  pre^ent^d  in  Chapter  7. 
The  actual  rocket  forces  and  pointing  angles  (f  9  f  attained  are 


determined  by  actuation 


on  the  idealized  commands 


system  transfer  functions  (G‘  ^  ,G^  ) 

.  a  .‘/S.  '  ' 


acting 


The  measures  for  evaluation  of  the  control  performance  can  be 
divided  into  two  distinct  classes:  (1)  the  performance  of  the  control 
law  in  steering  the  ejection  seat  along  a  prescribed  trajectory,  (2) 
the  performance  of  the  control  law  in  bounding  the  lethality  measure 
or  acceleration  radical  given  the  constraints  of  the  prescribed 
trajectory.  In  this  particular  application  the  "prescribed 
trajectory"  has  the  loose  meaning  of  ground  collision  avoidance, 
dynamic  pressure  reduction,  and  control  of  the  vehicle  to  terminal 
attitude  conditions  conducive  to  succeeding  phases  of  chute 
deployment.  The  acceleration  radical  employed  in  this  study  is  given 
by: 


( 


b  b 

X  v^ 


•ti  -b 

-  V.  +  wj  .  R  ) 


(6.1) 


The  radical  coefficients  (a)  in  equation  6.1  represent  the  relative 
sensitivity  to  injury  of  the  pilot  by  mutually  orthogonal  acceleration 
components  acting  normal  to  the  chest,  laterally  and  down  the  spine. 
The  specific  values  for  coefficients  (a)  expressed  in  the  square  of  g 


[(l/30g)^,  (1/I2g)^  (1/I7g  )^1 


The  acceleration  radical  is  assumed  to  be  lethal  or  at  least 
implies  a  high  probability  of  serious  injury  whenever  a^  >  1. 
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The  sinulaclon  segments  from  Figure  6. 1  mechanized  in  equivalent 
analog  circuitry  are:  (1)  Aerodynamic  Coefficients  (Block  Bl),  (2) 
6D0F  Flight  Model  (Block  B2),  (3)  Sensor  Models  (Block  B3),  (4) 
Actuator  Dynamics  Model  (Block  B6).  and  (5)  Resultant  Evaluation 
(Block  B7).  The  ensuing  description  (paragraphs  6.2. 1-6. 2. 6)  dlcusses 
Che  mathematical  models  employed  and  defines  Che  parameters  utilized. 
Digital  operations  is  reserved  for  mechanizations  of  (6)  Control 
Algorithm  (Block  B4)  and  (7)  Idealized  Rocket  Configuration  (Block  B5) 
discussed  in  paragraphs  6. 3. 1-6. 3. 2. 

6. 2  Analog  simulation  models. 

6.2.1  Aerodynamic  Model.  The  source  of  aerodynamical 
coefficients  for  the  ejection  seat  is  taken  from  the  report  entitled 
"Aeromechanical  Properties  of  Ejection  Seat  Escape  Systems",  B.J. 

White  (1974),  That  report  dealt  with  Introductory  technical 
discussions  on  aerodynamical  coefficients,  forces  and  moments  with  an 
emphasis  on  presenting  Che  measuTements  of  aerodynamical  coefficients 
for  crew  escape  systems  performed  at  the  Air  Force  Flight  Dynamics 
Laboratory,  Wright  Patterson  AFB.  Intensive  aerodynamic  measurements 
of  the  F-lOl  and  F-106  manned/ unmanned  ejection  seat  systems  were 
performed  over  a  wide  range  of  aerodynamic  conditions  (0.2-1, 5  Mach, 
0-360  deg  attack  angle,  0-45  deg  sideslip  angle)  resulting  in  the  most 
extensive  compendium  of  such  data  to  our  knowledge. 


The  B.J.  White  data  for  the  manned  ACES  ejection  seat  form  the 
basis  for  Che  aerodynamical  model  for  the  results  presented  herein. 

The  raw  aerodynamical  coefficients  are  defined  In  the  following  cable: 


'JV L^  L-V  i.-*.  1  ■'-  -W  '.-^ 


Table  6.2  TefinUion  of  fundamental  Aerodynamical  Toef  f  Ic  iencs. 
Coefficient  Definition 

axial  force  coefficient 

Cy  side  force  coefficient 

C  j  normal  force  coefficient 

C  rolling  moment  coefficient 

C  ^  pitching  moment  coefficient 

C  ^  yawing  moment  coefficient 

roll  acceleration  damping  coefficient 

C  pitch  acceleration  damping  coefficient 

ra 

C  yaw  acceleration  damping  coefficient 

"r 

The  aerodynamical  forces  (f^  )  and  torques  ( )  expressed  in 

the  ejection  seat  body  frame  are  related  to  the  reference  coefficients 
in  Table  6.2  by: 


-  Q  s 


t‘.  g. 


C^(a.3i 

cy(5,3,  n) 

L  Cz(i,3,  n)J 


C^(^3,  n)  >  p  Cip  ^  Cy  y^c 

(a.  3.  n)  >  q  Cn,q  ^CG  ^z~  ^CG  ^  x 


C_(a,3.  n)  +  r  dhZ  ^  y^c  C^- 


2V 


The  aerodynamic  pressure  Q  is  defined  by  the  familiar  expression: 


Q  ■  4  0  (h,T)V 
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where  o  aCmospherlc  density  at  altitude 

h  altitude  of  vehicle 
T  reference  temperature  at  altitude 

The  aCmospherlc  density  as  a  function  of  altitude  Is  generated  by 
look-up  tables  equivalent  to  the  method  employed  In  the  EASIEST 
simulation  model  (1980). 

Actual  generation  of  aerodynamical  coefficients  for  Intermidlary 
points  lying  between  those  given  by  White  are  generated  by  linear 
Interpolation  over  the  appropriate  Mach, a  ,6  interval  since 
justification  of  a  higher  order  model  is  not  apparent  from  inspection 
of  the  data.  The  relative  Insensitivity  of  the  B.J.  White 
coefficients  over  Che  0.5  -  1.2  Mach  regime  is  Indicative  of  turbulent 
aerodynamic  flow  for  most  flight  conditions  of  interest  in  this  study. 
Further,  it  should  be  noted  that  the  control  synthesis  method 
discussed  in  Volume  I  is  highly  reliant  on  the  sensor  Inputs  while 
avoiding  any  explicit  estimation  schemes  dependent  on  aerodynamical 
coefficients.  Hence,  aerodynamic  effects  are  created  as  disturbances 
in  the  control  algorithm  which  manifest  their  influence  as  undesirable 
state  dynamics  Co  be  damped  out  by  considering  Che  variations  in  the 
observables. 

6.2.2  Six  Degree  of  Freedom  Flight  Model  (6  DOF).  The  6  DOF 
flight  equations  are  the  differential  equations  of  motion  representing 
the  evolution  of  the  system  states  in  body  coordinates  according  to 
Newton's  laws  of  motion.  The  inertial  forces  and  torques  acting  on 
the  body  are  Identifiably  from  three  distinct  sources:  (1) 
aerodynamic  forces  and  torques  (paragraph  6.2.1),  (2)  resultant  rocket 
forces  and  torques  (paragraph  6.2.5)  and  (3)  gravitational  force.  All 
ocher  terras  contributing  to  the  translational  and  rotational  state 
dynamics  are  merely  consequences  of  the  particular  frame  selected  for 
Che  mechanization.  Selection  of  Che  body  frame  for  the  representation 
of  rotational  motion  significantly  simplifies  Che  differencial 


equations  since  in  chat  frame  the  moment  of  inertia  is  constant. 

Hence,  it  is  Che  natural  frame  of  choice  for  implementation  of  control 
strategies  and  modeling  of  the  dynamics  of  the  rigid  body.  Excellent 
presentations  of  the  basic  theory  of  dynamics  may  be  found  in 
Goldstein  (1981)  or  Eckln  (1972). 

The  non-linear  6D0F  equations  in  the  ejection  seat  body  frame  for 
this  problem  are  summarized  in  Figure  6.2.  The  symbols  in  Figure  6.2 
were  previously  defined  in  Table  6.1.  Notice  that  the  first  six 
equations  in  Figure  6.2  are  the  expression  of  Newton's  laws  in  the 
rotating  body  frame  while  equations  7-9  are  Che  inertial  velocity 
components  expressed  as  functions  of  the  body  velocity.  Finally, 
equations  10-12  are  the  kinematic  equations  for  the  inertial  Euler 
angles  as  a  function  of  body  angular  rates  and  orientation.  The 
essential  singularity  in  equations  10-12  occures  at  9  *90  degrees 

due  to  the  order  of  definition  of  the  Euler  angles  so  that  the  pitch 
angle  limit  in  examining  control  law  performance  is  Inherently  bounded 
to  less  chan  the  straight-up  direction. 

6.2.3  Sensor  ftodels.  The  "acceleration  control"  approach 
requires  chat  rapid  reliable  estimates  of  inertial  forces  and  torques 
acting  on  the  ejection  seat  be  available  for  generating  immediate 
control  terras  that  neutralize  undesirable  state  dynamics.  In  addition 
the  control  approach  is  reliant  upon  direct  estimates  of  the  inertial 
Euler  angles  for  generation  of  control  torques  proportional  to 
attitude  error  in  order  to  steer  Che  ejection  seat  to  a  desired 
terminal  attitude.  The  requirement  to  directly  estimate  inertial 
torque  Imposes  at  a  minimum  that  measurements  of  Inertial  rotational 
rate  and  acceleration  be  continuously  available.  Rate  gyros  and 
angular  accelerometers  are  standard  inertial  instruments  that  provide 
these  measurements  and  are  small  in  size,  low  in  power  with  large  mean 
time  between  failures.  Similar  considerations  apply  for  the  estimate 
of  inertial  forces  and  velocities  imposing  the  additional  requiremant 
of  three  axis  accelerometer  readout.  The  attitude  estimates  are 


Figure  6.2  Six  Degrees  of  Freedom  Differential  Equations 


either  directly  avalleble  as  synchro  resolver  outputs  oi  Jciivable 
from  the  integrated  outputs  of  the  rate  gyros  with  proper 
initialization  prior  to  ejection  time  from  navigational  estimates  from 
the  host  aircraft.  Volume  I,  section  5.1  summarizes  the  results  of  a 
literature  survey  of  currently  available  Inertial  Instruments  that 
fulfill  the  sensor  requirements.  The  objective  In  Volume  II  Is  to 
relate  the  mathematical  form  of  the  sensor  models  and  to  establish  the 
level  of  error  Implied  by  the  measurements. 

6.2.3.  1  Sensor  Atcuracy  Limitations.  The  hybrid  simulation 
environment  necessitates  the  use  of  interfacing  electronics  that 
samples  the  continuous  analog  representation  of  input  variables  and 
delivers  the  quantized  "snapshot"  estimate  to  the  digital  control 
algorithm.  The  analog  to  digital  converter  (ADC)  is  the  standard 
device  for  such  purposes.  The  ADC  accepts  as  input  an  analog  signal 
specified  over  a  limited  range  in  voltage  magnitude  and  outputs  the 
digital  signal  equivalent  in  the  form  of  the  number  of  quanta 
expressed  in  the  base  unit  of  the  device.  The  maximum  number  of 
quanta  is  given  by  2  where  N  is  the  number  of  bits  of  accuracy  in 
ADC  resolution  so  that  the  base  unit  is  1/2  '  if  the  magnitude  of  the 
signal  is  bounded  by  0  and  1. 

The  actual  accuracy  delivered  by  ADCs  is  dependent  on  the 
electrical  noise  environment  where  it  operates  and  may  in  fact  be 
substantially  less  than  the  N  bits  of  precision  guaranteed  due  to  the 
presence  of  stray  electrical  noise  at  the  analog  interface.  The  ADCs 
used  here  are  rated  as  12  bit  devices  while  actual  noise  measurements 
at  the  a-/-  10  V  Interface  indicate  ambient  electrical  noise  at  the  ^0 
mV  level.  Hence,  at  best  the  ADCs  can  be  expected  to  deliver  10  bits 
of  precision  for  a  base  unit  of  i/2  rather  than  the  1/2  *■  rating. 

Table  6. 3  specifies  each  input  to  the  control  algorithm 
illustrated  in  Figure  6.1  with  associated  maximum  input  analog  bounds, 
corresponding  real  scale  factor,  base  unit  assumed  and  resulting 


quantization  error 


Table  6.3  Inertial  Instrument  Accuracy  Limits. 


Sensor  Variable  Analog  Voltage  Real  Scale  Base  Unit  Quantization 


Bounds 

Factor 

Error 

Synchro 

Yaw 

-►/-  10  V 

*!-  200  deg 

1/1024 

*l~  0.2  deg 

Resolver 

Synchro 

Pitch 

-*■/-  10  V 

•►/-  200  deg 

1/1024 

+/“  0.2  deg 

Resolver 

Synchro 

Roll 

■*•/-  10  V 

*l~  1000  deg 

1/1024 

+/-  1.0  deg 

Resolver 

Accelerometer 

+/-  10  V 

+/“  2000  fps^ 

1/1024 

+/-  2  fps^ 

Rate  Gyro 

+/-  10  V 

->•/-  10  R/s 

1/1024 

+/-  0. 5  deg/s 

Angular 

+/-  10  V 

2 

*!-  500  R/s 

1/1024 

2 

■♦•/-  30  deg/s 

Accelerometer 


It  Is  evident  from  Table  6.3  that  the  quantization  error  alone  would 
qualify  the  sensors  as  crude  when  compared  with  the  high  accuracy  of 
modern  inertial  instruments. 

A  more  complete  discussion  of  the  hardware  Interface 
specifications  for  the  hybrid  simulation  is  presented  in  paragraph 
6. A.  This  segment  is  included  to  quantify  the  sources  of  error  beyond 
the  Intended  sources  of  error  explicitly  modeled. 


6. 2.3.2  Sensor  Dynamics  Model.  For  inertial  navigational 
applications  It  Is  usual  to  consider  sources  of  error  in  sensor 
modeling  which  will  cause  apparent  velocity  and  position  error  growth 
when  integrated  over  the  long  term  (I  hour).  For  accelerometers  and 
gyros  the  major  sources  of  error  include  constant  bias,  random  bias, 
scale  factor  error,  non-linearity  of  the  scale  factors  and 
misalignment.  Fnr  typical  high  accuracy  inertial  instruments  in  a  low 
dynamic  environment  the  resultant  navigation  errors  are  generally  in 
the  vicinity  of  1  nm  position  error  and  2  fps  velocity  error  after  an 
hour  of  operation. 

Given  the  extremely  brief  interval  of  this  application  (at  most  2 
sec)  many  of  the  mentioned  sources  of  instrument  error  are  negligible. 
In  light  of  the  discussion  of  the  previous  paragraph  the  sources  of 
error  associated  with  bias,  g-sensitivlty ,  scale  factors  and 
misalignment  are  moot  when  compared  with  the  quantization  error  in 
Table  6. 3. 

Of  far  greater  consequence  in  the  ejection  seat  problem  is  the 
delay  Imposed  by  the  sensors  given  the  reliance  of  the  control 
strategy  on  the  availability  of  vehicle  state  estimates  and  the  need 
to  respond  immediately.  The  emphasis  on  sensor  modeling  in  this  study 
then  is  to  consider  the  effects  of  sensor  delay  on  the  control  system. 
The  highest  fidelity  models  for  sensors  generally  are  high  order  (up 
to  sixth)  which  generate  with  considerable  resolution  the  sensor 
dynamics.  In  the  interest  of  simplicity  the  delay  models 
utilized  here  are  of  first  order  which  is  meant  to  capture  the  essence 
of  system  delay  on  control  performance. 


Table  6.4  defines  the  relationships  between  the  sensor  estimates 

(  1  ,a^  ,v^  ,  d,^)  to  the  truth  values  ( i  ,a^  ,v^  ,  w*?  ,  )  as  well 

i-iii  ii-i-i- 

as  the  transfer  functions  (G.^  ,G  ,G  ,G-  )  with  the  associated 


time  constants  (  ,  ^a  ,  )•  Notice  that  all  transfer  functions 

I  1  i  i 

are  defined  in  the  frequency  domain  while  the  error  in  the  velocity 
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estimate  v^'^s)  Is  derived  from  the  acceleration  error  (  a^(s)  - 

3*^(8)).  The  time  constants  In  Table  6.4  are  representative  values  for 
1 

sensors  deployed  In  missile  guidance  applications.  It  should  be  noted 
that  the  control  scheme  was  evaluated  with  the  time  constants 
^3^,..,  sec  (well  within  current  sensor  bandwldths)  and  yielded 

results  similar  to  those  presented  in  Chapter  7. 

Table  6.4  Sensor  Transfer  Functions  and  Response  Times. 


1^(3)  -  Gj 

(s) 

a^  (s)  -  ^(s) 

v^  (s)  -  v^  (s) 

(s)  -  G^ 

(s)  -  g/ 
i 

(s) 

(s) 

b  y  V 

.  (s)  -  a 
1. 

\  (s)) 

Sensor 

Transfer 

Function 

Time  Constant 

Synchro  Resolver 

C 

1 

^3^-  0.00025 

Accelerometer 

^1 

G 

1  +  S 

1 

T  -  0.0025 

Rate  Gyro 

C(j)  ■ 

X 

1 

T  -  0.0025 

1  -t-  T  S 

“  i 

Angular  Accelerometer 

^  ” 

1 

1 

T.  -  0.0025 

-J  . 

1  +  TV  S 

1 

6.2.4  Actuator  Dynamics  Model.  The  actuation  system  elements  In 
ejection  seat  control  are  rocket  nozzle  thrusters  either  fixed  In 
orientation  with  respect  to  the  ejection  seat  or  glmballed  to  allow 
for  some  degree  of  freedom  In  force  vectoring.  Emerging  technologies 
In  rocket  propulsion  systems  also  allow  for  generation  of  variable 
force  magnitudes  from  the  rocket  elements.  The  summary  of  a 
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technology  survey  for  currently  available  actuators  was  presented  in 
Volume  I,  paragraph  5.5.2,  This  paragraph  Is  concerned  with  the 
response  model  of  these  rocket  elements  to  input  commands.  It  is 
apparent  that  the  actuation  system  elements  have  finite  bandwidth  and 
representing  the  delays  imposed  are  important  In  evaluating  the 
control  system  robustness  to  such  restrictions. 


The  idealized  actuator  configuration  discussed  in  detail  in 

Volume  I  consisted  of  three  rocket  thrusters,  all  with  with 

variable  force  iiagnitudes,  one  fixed  in  orientation  with  respect  to 

the  seat  (main  thruster)  and  the  ocher  two  gimballed  with  two  degrees 

of  freedom.  The  independent  command  input  variables  with 

configuration  total  seven:  three  force  magnitudes  (f.' f.  ,f-  ^  and 
. M2  S,.  .  ^ 


four  rocket  pointing  angles  i.e.,  two  angles  each 

C  C 

for^  two  rockets.  The  relationships  for  generating 


control  forces  (f®  ) 


,f^  -',■5^  ^  Idealized 

Q 

and  torques  (t,,)  indicated  in  Figure  6.1  are  developed  in  Volume  I, 
paragraph  3.2  and  summarized  in  paragraph  6.3.2  o:  this  volume. 


In  the  interest  of  ease  in  implementation  first  order  models  are 
incorporated  to  corrupt  the  idealized  input  commands  generated  by  the 
control  algorithm.  Table  6.5  relates  outputs  to  inputs  and  defines 
the  transfer  functions  and  response  times  utilized  in  generating  the 
achieved  actuation. 


Table  6.5  Actuator  Transfer  Functions  and  Response  Times. 


Actuator  Update 


S 


1 


e 


Nl 

9c 


Actuator  Transfer  Time  Constants 

Element  Function 


Pointing  Angle 


.Nl 


1-f  T 


eNi 


.Nl 


■  0.015  sec 


Force  Magnitude 


TTT - ;  i  *  0.0025  sec 


The  connanded  values  In  Table  6.5  are  quantized  consistent  with  the  12 
bit  DACa  and  real  scale  factors.  In  actuality  the  DAC  performance  Is 
limited  by  ambient  electrical  noise  at  the  output  channels  analogous  to 
noise  Input  corruption  for  the  sensors.  Direct  measurements  of  the  noise 
levels  Implied  at  best  11  bit  delivery.  Table  6.6  defines  the  output 
variables  along  with  real  scale  factor,  base  unit  and  the  resulting  maximum 
quantization  error. 
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Table  6.6  Actuator  Element  Accuracy. 


Variable 

Definition 

Rase  Unit 

Scale 

Factor 

Quantization 

Error 

CD 

0 

pitch  pointing  angle 

for  rocket  N’. 

I 

1/20A8 

+  XR 

+  0.2  degree 

yaw  pointing  angle 
for  rocket  Nj 

1/2048 

+  xR 

+  0.2  degree 

force  magnitude  for 
rocket 

1/2048 

+  lOOOO 

lb  +  10  lb 

6.2.5  Resultant  Evaluation.  The  6  DOF  flight  model  discussed  In 
section  6.2.2  required  the  resultant  inertial  forces  and  torques 
acting  on  the  ejection  seat  In  order  to  propagate  the  state  according 

to  the  differential  equations  of  motion.  The  Inertial  roclcet  force 
b  b 

(f^  )  and  torque  )  resultants  expressed  In  body  appear  directly  In 

the  6D0F  equations  as  drivers  of  the  system  states.  The  resultants 
b  b 

(fj.  ,'r  )  are  simply  the  summation  of  the  Individual  rocket  inertial 
components  transformed  to  body  coordinates.  As  sucn  they  may  be 
expressed  succintly  as : 


^  (0 


^'l 


N, 


5r 


X  C 


N. 

1 


-^1 

-  )f  (6.2) 


N  N 

N'.  1  -  i 

(e  )f  (6-3) 


In  equations  6. 2-6. 3  f 

b 


1  i 

,'3  ,  ‘  are  of  course  the  achieved 


actuation  signals  while  Cm  are  the  vectors  of  direction  cosines  from 

•  1  r  ^'i 

rocket  to  body.  The  parameter  Is  the  rocket  displacement 

from  the  seat  c.g.  expressed  in  body  coordinates.  Rockets  Ni .No  are 

glmballed  with  pointing  angles  (9  1,9  ‘^)  while  N3  Is  fixed 

b 

leading  to  expressions  for 


■CeNiC/i  ■ 

V 

1 

C  - 

1-1.2  ; 

L  tan  9  J 

Ls.N  C,  N, 

The  ootlvatlon  for  selection  of  this  particular  actuation 
configuration  and  the  development  of  the  theory  may  be  found  In  Volume 
I,  section  3.2. 

6.2.6  Summary  of  Analog  Simulation  Models.  The  preceding 
paragraphs  provided  a  synopsis  of  the  mathematical  form  of  the  analog 
models  while  omitting  completely  the  analog  schematics  which  Implement 
the  circuit  representation  of  those  models.  The  detailed  description 
of  the  analog  hardware  Is  contained  In  Volume  IV,  Appendix  F,  entitled 
"ACES  II  Hybrid  Simulation  Report"  prepared  by  the  Technical  Computing 
Center  of  MMOA.  Some  of  the  Information  presented  herein  are  excerpts 
from  the  material  contained  in  Appendix  F.  As  a  matter  of  correlating 
Che  Information  In  Volumes  II  and  IV  Table  6.7  defines  Che 


cross-reference  Index  relating  the  analog 

schematic  representations  o: 

Volume  IV  to 

Che  machematlcal  form  of  Volume  II. 

Table  6. 7 

Correspondence  of  Mathematical 

Model  to  Analog  Schematics 

Volume  II  Subject 

Volume  IV 

Paragraph 

Analog  Schematic  Form 

6.2.  1 

Aerodynamic  Model 

NDE 

6.2.2 

6D0F  Flight  Model 

pp.  115-121 

6.2.3 

Sensor  Models 

p.  121,124 

6.2.4 

Actuator  Dynamics  Model 

p.  122 

6.2.5 

Resultant  Evaluation 

p.  123,125 

NDE-no  direct  equivalent 


It  should  b«  noted  that  the  Aerodynamical  Model  In  Table  6.7  has 
no  direct  analog* equivalent.  The  evaluation  of  aerodynamic 
coefficients  Is  implemented  on  a  parallel  processor  (AD-IO)  configured 
with  an  (ADC/DAC)  interface  to  the  6D0F  flight  model. 

This  concludes  Che  discussion  of  the  Analog  Simulation 
requirements  for  the  Hybrid  Simulation  development.  The  digital  ^ 
segment  is  concerned  with  generation  of  Che  control  signals  it.  ) 

♦  A,  k  * 

given  the  sensor  Inputs  .a^^  .v^^  )  which  is  the  subject  of 

Che  next  paragraph. 

6. 3  Digital  Algorithms. 

The  digital  segment  of  the  hybrid  simulation  is  the  realization 
of  real-time  control  algorithms  which  process  sensor  measurements  and 
produce  control  commands  for  use  by  the  actuation  system  in  steering 
the  ejection  seat  along  a  prescribed  trajectory.  The  rapid  update 
race  (SOHz)  mandates  the  use  of  numerically  efflcier'.c  algorithms  which 
evaluate  inertial  torques  and  forces  and  apply  corrective  terms  In 
order  to  meet  Che  trajectory  specifications. 

The  sensor  l-oucs  were  dealt  with  In  paragraph  6.2.3  and 
paragraph  6.2.4  defined  Che  dynamics  of  Che  actuation  system  elements. 

The  objective  here  Is  to  summarize  Che  processing  by  the  control 

b 

algorithms  that  generate  Che  Idealized  body  forces  )  3^(1  torques 
(  )  and  transform  them  Into  desired  actuation  signals. 


6.3.1  Control  Algorithm.  Development  of  Che  control  system 
methodology  was  a  major  topic  of  Volume  I.  In  chat  discussion  the 
exigency  of  Che  ejection  seat  control  problem  was  emphasized  which  led 
CO  Che  "acceleration  control"  concept.  The  conclusion  of  chat 
investigation  was  Che  necessity  for  prompt  neutralization  of  the 
translation  and  rotational  accelerations  which  led  to  specific  sensor 
and  actuator  system  requirements.  Given  that  the  aerodynamic  effects 


are  reduced  to  acceptable  levels  the  problem  of  meeting  terminal 
attitude  constraints  Is  a  residual  dynamics  problem  amenable  to 
standard  attitude  control  methods  such  as  "cr  ss  product  steering". 

The  translational  problem  is  best  dealt  with  by  selecting  desired 
acceleration  profiles  and  generating  the  associated  steering  commands 
to  meet  the  profile  constraints.  The  final  control  logic  resulting 
from  the  above  considerations  is  detailed  in  Volume  I,  paragraph  4.3. 
In  addition  the  mechanization  of  that  control  logic  is  also  summarized 
in  Volume  IV,  pp.  8-13. 

The  control  logic  is  briefly  reviewed  here  to  allow  for 
continuity  in  the  discussion.  Given  the  terminal  attitude  constraints 
(  )i  Che  intent  is  to  steer  the  ejection  seat  from  some  attitude 

(  $  (  9  )  to  the  terminal  conditions.  Small  angle  errors  may  be 
corrected  effectively  by  "cross  product  steering".  In  that  method 
torque  terms  proportional  to  angular  error  and  its  rate  are  returned 
as  feedback  to  achieve  and  maintain  the  desired  regulated  attitude. 

Figure  6.3  sketches  the  attitude  control  logic  for  cross  product 

steering.  The  objective  is  to  align  the  unit  vertical  body  vector  z^ 

with  the  terminal  desired  body  vector  zWji  9J.The  error  signal  is 
b  b  F  F  f  r 

generated  by  Zg  X  ,  exact  for  small  attitude  deviations  and  hence 
the  use  of  the  term  "cross  product".  Step  3  in  Figure  6.3  defines  the 
evaluation  of  the  desired  angular  acceleration  proportional  to 
attitude  error  and  its  rate.  That  is  kgpg  is  the  position  gain  while 
k  CPj^  serves  as  the  velocity  gain.  Step  4  determines  the  desired 
rotation  rate  from  the  angular  acceleration  scaled  by  "time  to  go"  as 
an  approximation  to  the  exact  integral.  The  nonlinear  form  of  the 
control  is  evident  in  step  5  which  also  directly  cancels  the  inertial 
torques  acting  on  the  seat  ^^X  I^^  ).  The  expression  (5)  is 

the  major  driver  for  attitude  steering  while  the  linear  contributions 
in  (6)  become  apparent  for  large  deviations  from  the  nominal 
trajectory. 


Step  I : 


Step  2: 


Step  3: 


Step  4 ; 
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Step  5 ; 
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Step  6 : 


1.  Inertial  to  body  direction 
cosine  evaluation 


2.  Calculation  of  final 

attitude  vector  In  body 


3.  Evaluation  of  desired 
angular  acceleration  In 
body 


4.  Projection  of  desired 
angular  rate  In  body 
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6.  Delta  Torque  correction 
feedback. 


7.  Final  Inertial  torque  In 
body  coordinates 


Figure  6.3  Attitude  Steering  Logic  for  Cross-Product  Law 


The  brief  description  above  related  to  the  small  angle  error  case 
while  the  considerations  for  the  large  angle  case  is  similar  with  some 
modifications  in  steps  2-3.  The  large  angle  error  case  by  comparison 
considers  only  the  angular  error  in  step  3  to  arrive  at  the  desired 
angular  acceleration  while  preserving  the  form  (4-7),  The  nonlinear 
compensation  (5)  performs  similar  to  linear  control  for  small 
deviations  from  the  nominal  trajectory  while  generating  quadratic 
correction  feedback  of  acceleration  for  angular  rate  error  in  response 
to  linear  growth. 

The  translational  logic  is  illustrated  in  Figure  6.4,  The 

desired  inertial  velocity  v  ^  is  analytically  the  integral  of  the 

*^1  i  i 

desired  inertial  acceleration  a^  .  The  dynamic  values  for  v^  ,  a^  in 

turn  reflect  the  initial  and  terminal  desired  acceleration  and 

velocity  conditions  to  which  they  converge  over  the  ejection  control 

Interval.  The  variables  v^  ,  a^  at  ejection  time  are  the  initial 

d  d 

values  of  velocity  and  acceleration  and  are  selected  to  change  slowly 
in  time  approaching  the  terminal  conditions.  As  a  result  a  smooth 
transition  from  the  initial  to  the  final  desired  values  is  realized. 

The  commanded  forces  are  generated  in  the  body  frame  while  the 

updates  of  a ^  and  Vj  are  strictly  inertial  (step  2).  Expression  (3) 
d  d 

cancels  acting  inertial  forces  (-m  a.  )  and  supplants  the  same  with 

^  b 

the  desired  dynamics  imbedded  within  the  commanded  acceleration  a^  • 

Expression  (4)  is  proportional  and  integral  compensation  for  error 

drift  from  the  nominal  velocity  and  acceleration  profile  resulting  in 

the  dominance  of  a  j  in  the  terminal  conditions, 
d 

n  g 

The  feedback  forces  (f  )  and  torques  (t  )  of  Figures  6. 3-6. 4 

c  c 

form  the  idealized  commands  in  the  body  frame  and  are  synonymous  with 
those  defined  in  Figure  6.1.  In  any  implementation  of  controllers  for 
ejection  seats  the  idealized  commands  need  be  converted  to  actuation 
system  element  commands.  The  discussion  of  the  actuation  system 
mechanization  is  the  subject  of  the  next  paragraph. 


1.  Desired  Velocity,  force 
In  body  from  Inertial 
components 


Step  2. 


2.  Update  of  desired  inertial 
acceleration  and  velocity 
components 


Step  3  : 


3.  difference  In  inertial 
rorce  and  desired  force 


Step  4  : 


4.  Delta  Force  Correction 
Feedback 


Step  5 : 


5.  Final  Inertial  force  In 
body  coordinates 
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Figuro  6  4  Force  Steering  I. ogle  in  Translational  Control 


6.3.2  Idealized  ActuaCor  Conf Iguraclon  Model.  The  actuation 
system  employed  is  discussed  in  detail  In  Volume  I,  section  3. 2. 2. 3. 
Similarly  the  processing  steps  are  summarized  in  Volume  IV,  pp. 

13-14.  The  excerpt  from  Volume  IV  is  included  to  reduce  the  necessity 
of  cross-reference  in  pursuing  this  discussion. 


D  Q 

The  problem  here  is  to  convert  the  idealized  commands  (fg  ,t  ) 

—  Nj  m,  ^ 

into  the  nearest  equivalent  set  of  rocket  commands  (f  ,0  ^). 

N  c  ’  c 

term  T  is  the  rocket  force  magnitude  while  9  , i  are  the  rocket 

c  c  c 


The 


pitch  and  yaw  pointing  angles  respectively.  Previously  defined  in 


paragraph  6. 2. 3. 2,  there  are  a  total  of  seven  control  elements  to  be 
_Nl  -^2  -^3  S.  No  No 

considered:  (f  ,f  ,f  ,0  ,  i.e.  ,  three  force 

^  c  c  c 


magnitudes  and  four  rocket  pointing  angles. 


Conceptually  the  problem  may  be  seen  as  choosing  the  parameters 

_Ni  Ni  Si  . 


(fj,  ,0  ,vc  )  to  simultaneously  meet  the  decoupled  idealized  commands 


(f®  ,0)  and  (0,T®  ).  That  is,  rocket  commands  allocated  to  meet  the 


resultant  force  f®  will  generate  no  disturbance  torques.  An 


analogous  statement  applies  for  the  command  pair  (0,T®  )  so  that  the 

c 


rockets  would  deliver  exact  force  command  pair  cancellation  resulting 
in  no  effect  on  CG  translational  motion. 


Figure  6.5  defines  the  relationships  for  producing  the  six 

dimensional  force  vector  (f,,...,f  )  from  any  given  command  pair 

^  o 

(f?  I T®  )•  The  terms  Q.  ,  are  the  inverse  terms  of  the  control 
c  c  ij 

gradient  matrix  P  which  relates  each  rocket  contribution  to  the 

command  pair  resultant.  Once  available  the  force  vector  (f,  ,...,f,  ) 

-N.  N.  i  0 

are  converted  to  the  actuator  signals  (f_^,9  ,  ;  as  defined  in 

c  c  c 

Figure  6.6.  There  is  an  apparent  discrepancy  in  dimension  between  the 
input  vector  (dimension  six)  and  the  output  vector  (dimension  7).  The 
extra  degree  of  freedom  is  removed  by  the  arbitrary  assignment  of 
rocket  1,2  sharing  of  element  f^  evident  in  expressions  1-2  of  Figure 
6.6.  This  assignment  also  affects  the  generation  of  the  rocket 
pointing  angles  in  expressions  (5,7). 


>  inf  ill}* 


I'aluation  of  Idealized  Command 


The  generation  of  the  rocket  commands  completes  the  description 
of  the  processing  requirements  by  the  digital  algorithms.  The 
hardware  allocation  to  support  the  processing  of  the  Analog  and 
Digital  models  is  described  next. 

6. 4  Hardware  description. 

6.4.1  Requirements.  At  a  minimum  the  hybrid  simulation  requires 
the  use  of  analog  computing  facilities  for  implementation  of  the 
analog  models  of  section  6.2,  the  development  of  a  microprocessor 
based  hardware  controller  with  associated  Interfacing  electronics  plus 
some  means  for  displaying  the  state  variables  and  system  performance 
measures.  The  hybrid  computing  facility  at  MMOA  has  at  its  disposal  a 
powerful  array  of  parallel  processors,  digital  and  analog  computers 
and  chart  recorders  for  hosting  complex  hybrid  simulations. 

Linidynamics  of  Phoenix  Arizona  (UPHX)  was  tasked  with  the  development 
of  the  hardware  controller  board  that  supported  the  interface  and 
computational  requirements  for  the  Motorola  68000  microprocessor.  In 
this  phase  of  the  effort  SSI  provided  the  system  specification  and  was 
responsible  for  development  of  the  microprocessor  based  control 
algorithms. 

Figure  6.7  illustrates  the  division  of  responsibility  for  the 
implementation  of  the  integrated  simulation  and  control  system. 
Inspection  of  Figure  6.7  indicates  the  hardware  allocation  to  meet 
system  requirements.  The  following  paragraphs  describe  the  various 
processors  in  use. 

6.4.2  MMOA  Analog  Hardware  Description.  Volume  IV  details  the 


hardware  sub-systems  and  applications  software  necessary  for  operation 
of  the  analog  segment  of  the  system  in  Figure  6.7.  A  brief  synopsis 
of  those  hardware  elements  is  given  here  for  completeness. 


Table  6.8  Definition  of  Processors  in  MMOA  Analog  Simulation. 


Hardware  Definition  Quantity 
Component 

AD-10  Applied  Dynamics  AD-IO  Parallel  Processor  1 
EA-8800  Electronic  Associates  8800  Analog  Computer  3 
PE  8/32  Perkin  Elmer  8/32  Digital  Computer  1 
Strip  Chart  Recorders  7 


The  three  EA-8800  analog  computers  support  the  implementation  of  all 
analog  models  with  appropriate  ADC/DAC  Interfaces  for  communication 
with  the  AD-10,  PE  8/32  and  the  Unidynamics  controller  board.  The  PE 
8/32  defined  in  Table  6.8  does  not  explicitly  appear  in  Figure  6.7. 

The  MMOA  analog  simulation  allows  the  use  of  the  PE  8/32  to  serve  as 
substitute  in  place  of  the  hardware  controller  microprocessor  based 
system  under  test.  Section  6.4.4,  Modes  of  Operation,  describes  the 
use  of  this  alternate  system  for  the  generation  or  results.  The  strip 
chart  recorders  allow  for  visibility  of  dynamic  v-iriables  each  with 
preset  maximum  values  divisible  over  a  100  to  1  scaling  range 
facilitating  the  observation  of  small  signal  dynamics  for  key 
variables  of  interest. 

6.4.3  Unidynamics  Hardware  Controller  Board  Description.  The 
general  requirements  for  the  fabrication  of  the  microprocessor 
controller  board  are:  (1)  Interface  electronics  to  allow  for  access 
of  sensor  Inputs,  (2)  sufficient  RAM  memory  to  support  the  control 
algorithms,  (3)  a  system  clock  with  a  suitable  rating  to  support 
real-time  operation,  (4)  Interface  electronics  to  allow  for  control  of 
actuation  system  elements  and  (5)  a  communications  Interface  to 
support  file  downloads  to  the  68000  microprocessor. 


Figure  6.8  lllustraces  che  syscem  specification  for  the 
Unldynaoics  hardware  controller  board.  The  system  Input  module 
consists  of  a  parallel  16  channel  12  bit  ADC  Interface.  The 
channels  receive  the  raw  analog  signals  from  the  MMOA  simulation  and 
deliver  che  digital  equivalent  to  Che  algorithms  hosted  on  che 
microprocessor.  The  system  outputs  are  serviced  by  a  parallel  12 
channel  12  bit  DAC  interface  with  sample  and  hold  circuits.  The 
assigned  control  syscem  variables  are  indicated  In  Figure  6.8.  Note 
Chat  che  current  configuration  requires  7  channels  to  service  outputs 
allowing  for  expansion  in  modeling  alternate  actuation  system 
configurations. 

The  RS-232  terminal  port  allows  for  direct  comaunlcatlon  with  the 
development  syscem  (VAX)  to  support  fast  download  of  executable  files 
in  che  prescribed  Motorola  format.  The  internal  memory  allocation  of 
the  microprocessor  allows  for  volatile  all  RAM  operation  (64K)  or 
permanent  program  storage  with  che  non  volatile  ROM  (64K)  memory. 

The  operation  of  the  sub-system  given  in  Figure  6.8  proceeds  as 
follows.  Ac  che  stare  of  an  execution  cycle  each  sensor  input  is 
requested  sequentially  over  the  parallel  input  communication 
interface.  With  appropriate  scaling  of  che  sensor  inputs  the 
controller  performs  calculations  to  generate  che  control  commands  for 
che  Individual  control  elements.  Each  of  the  digital  controls  to  the 
actuation  syscem  elements  are  converted  sequentially  by  che  parallel 
output  interface  to  analog  form  for  use  by  subsequent  analog  models. 
Hence  che  inputs  capture  all  state  estimates  at  nearly  che  same 
instant  in  time  and  che  outputs  are  available  for  che  actuation  syscem 
as  soon  as  practicable  with  near  simultaneity.  The  sample  and  hold 
circuitry  maintains  che  control  signals  at  che  commanded  level  until 
refreshed  by  che  ensuing  control  cycle. 
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The  RS-232  terminal  to  microprocessor  communication  Interface 
allows  for  the  control  algorithm  to  be  interrupted  In  order  to  view 
the  status  of  input  and  output  variables  or  internal  dynamic 
quantities  in  debugging  when  problems  do  arise. 

6.4.4  Modes  of  Hybrid  Simulation  Operation.  The  MMOA  analog 
system  described  above  alluded  to  direct  use  of  a  PE  8/32  digital 
computer  to  support  control  system  functions  as  an  alternative  to 
microprocessor  hardware  In  the  loop  operation.  The  PE  8/32  operating 
mode  executes  an  all  FORTRAN  Implementation  of  the  control  concept 
replacing  the  digital  control  segment  In  Figure  6.7.  This 
mechanization  Is  extremely  useful  as  It  allows  for  direct  comparisons 
with  microprocessor  In  the  loop  results  and  simplifies  alterations  in 
the  control  logic  for  the  evaluation  of  alternative  control  designs. 

In  addition  the  PE  operation  mode  allows  for  advanced  debugging  for 
isolation  of  potential  sources  of  system  mal-f unctions ,  both  hardware 
and  software. 

The  second  mode  of  hybrid  operation  is  of  course  with  the 
microprocessor  controller  in  the  loop.  This  mode  of  operation  Is 
primarily  concerned  with  verification  of  proper  operation  of  the 
microprocessor  based  code  and  to  provide  for  throughput  estimates. 
While  the  variables  assigned  to  the  ADC/DAC  Interface  are  directly 
observable  in  this  operating  mode  the  internal  controller  dynamic 
variables  are  unobservable  due  to  limited  communication  so  that  the  PE 
8/32  does  provide  for  a  more  complete  picture  of  controller  operation. 
With  that  in  mind  the  stand  alone  PE  8/32  operating  mode  is  actually 
preferable  for  the  examination  of  the  results. 


CHAPTER  7 
RESULTS  AND  ANALYSIS 


7 . 1  Overview  of  Chapter  7, 

Chapter  7  examiaes  the  results  obtained  with  the  proposed  control 
algorithm  when  exercised  in  the  real  time  hybrid  simulation 
environment.  The  major  goal  in  all  cases  is  to  demonstrate  acceptable 
performance  of  the  control  concept  with  respect  to  the  following 
measures:  (1)  avoidance  of  ground  collisions  for  low  altitude 

ejection,  (2)  containment  of  the  acceleration  radical  to  survivable 
limits  and  (3)  guarantee  that  in  steady  state  that  the  terminal 
attitude  constraints  are  met.  The  cases  under  investigation  vary  over 
a  wide  range  of  the  flight  envelope  expected  in  emergency  ejection. 
These  conditions  include  low  altitude  ejection  at  high  dynamic 
pressure,  low  altitude  ejection  at  "adverse  attitude"  and  high 
altitude,  high  dynamic  pressure  ejection.  It  is  important  to 
emphasize  that  the  results  discussed  below  utilize  a  single  set  of 
control  parameters  in  the  evaluation  of  all  the  test  conditions.  That 
is,  no  specific  tuning  of  parameters  to  the  specific  ejection 
conditions  are  inherent  in  the  following  results.  Any  viable  scheme 
necessarily  would  only  require  a  minimal  set  of  parameters  to  be 
selected  if  the  approach  is  truly  applicable  to  a  wide  region  of  the 
ejection  envelope. 

7.2  Comparison  of  Alternative  Operational  Modes  of  the 

Hybrid  Simulation. 

Three  operational  modes  are  examined  with  the  hybrid  simulation 
model.  The  first  illustrates  the  negative  consequences  of  open  loop 
ejection  where  there  is  no  direct  control  of  the  life  threatening 


forces  and  torques  In  ejection  nor  are  any  position  or  attitude 
constraints  enforced.  In  this  simplistic  mode  of  operation,  the  main 
thruster  is  deployed  at  a  constant  thrust  level  throughout  the 
scenario  and  all  other  control  elements  are  constrained  off.  The 
second  mode  of  operation  exercises  the  control  logic  hosted  on  a 
mainframe  computer  (PE  8/32)  which  serves  as  substitute  for  the 
microprocessor  based  controller  (see  paragraph  6.4.4).  In  this  mode 
of  operation  all  control  functions  are  exercised  to  simultaneously 
contain  the  acceleration  radical  and  to  provide  corrective  attitude 
steering  to  the  terminal  constraint  values.  The  results  obtained  with 
the  PE  8/32  represent  the  nominal  performance  to  be  expected  when  the 
controller  is  deployed  actively  in  emergency  ejection.  In  the  PE  8/32 
mode  of  hybrid  operation  the  control  logic  is  implemented  in  the 
FORTRAN  language  and  is  a  wholly  real  time  operational  mode.  The 
third  mode  of  hybrid  operation  is  with  the  control  law  hosted  on  a 
Motorola  68000  microprocessor  which  exercises  all  control  functions 
and  represents  the  performance  to  be  expected  with  actual  hardware  in 
the  loop  operation  with  accuracy  limited  to  16  bit  precision. 

In  summary  it  would  be  expected  that  the  open  loop  mode  of 
operation  would  typically  fall  with  respect  to  all  performance 
measures  for  "difficult"  ejection  situations.  The  PE  8/32  and  68000 
implementations  are  expected  to  deliver  acceptable  performance  with 
nearly  identical  results  since  both  modes  of  operation  are 
realizations  of  hardware  implementations  of  the  control  concept  under 
test . 


7.2.1  Interpretation  of  Strip  Chart  Results. 

The  outputs  of  the  hybrid  simulation  are  most  easily  assimilated 
as  direct  strip  chart  outputs  of  the  internal  analog  variable  values. 
A  typical  performance  evaluation  strip  chart  is  illustrated  in  Figure 
7.1.  The  key  variables  of  interest  in  evaluation  of  the  performance 
measures  are  included  in  the  figure.  Other  output  variables  are 


available  from  the  hybrid  simulation  but  are  omitted  in  the  interest 
of  summarizing  performance  in  limited  space.  Each  strip  chart  is 
divided  into  three  sequential  evaluations,  one  for  each  mode  of 
operation  identified  in  the  last  paragraph.  In  Figure  7.1,  each 
execution  mode  is  identified  at  the  top  of  the  figure  by  the 
abbreviations  "U","S",  and  "H".  Case  "U"  is  the  "uncontrolled" 
response  of  the  ejection  seat  and  represents  the  performance  of  fixed 
thrust  schemes  with  the  omission  of  STAPAC  from  the  ACES  II  ejection 
seat.  Case  "S"  outputs  represents  "simulated"  controller  performance 
and  are  the  results  of  the  PE  8/32  mode  of  operation.  Finally  the 
case  "H"  results  are  the  outputs  obtained  with  the  Motorola  68000 
"hardware"  in  the  loop  mode  of  operation.  All  of  the  results  examined 
in  this  section  are  of  2  seconds  duration. 

For  the  particular  results  illustrated  in  Figure  7.1,  the  initial 
flight  conditions  are  as  follows:  (1)  initial  speed  =  687  KEAS,  (2) 
roll, pitch, yaw  =  (0,0,0),  (3)  euler  rates  p,q,r  =(0,0,0)  at  (4) 
altitude  100  feet.  The  desired  terminal  conditions  desired  are  given 
by:  roll , pitch, yaw=(0,  10,0  deg)  with  ground  collision  to  be  avoided. 
The  ambient  dynamic  pressure  in  this  case  is  approximately  1600 
lb/ft**2  representing  a  high  dynamic  pressure  ejection  resulting  in  a 
lethal  condition  if  the  resultant  inertial  forces  and  torques  are 
inadequately  controlled. 

In  general  the  roll  euler  angle  is  essentially  a  "free"  control 
variable  since  rotations  about  the  velocity  vector  induce  no  injurious 
force  and  torque  contributions  of  any  consequence.  The  pitch  euler 
angle  is  subject  to  enormous  environmental  torquing  dynamics  so  that 
it  is  imperative  that  tight  pitch  control  be  maintained  to  Insure 
ejection  seat  body  stability.  Particularly  threatening  are  changes  in 
the  yaw  euler  angle  from  the  initial  condition  which  cause  large 
lateral  force  components  to  act  directly  along  the  most  sensitive 
pilot  body  direction.  It  is  apparent  that  the  yaw  euler  angle  should 
be  tightly  regulated  to  the  initial  value  in  order  to  minimize  the 


growth  of  the  lateral  acceleration  component. 

Examination  of  the  results  for  the  uncontrolled  ejection  seat 
labelled  "U"  in  Figure  7.  1  indicates  no  hope  of  pilot  survival  for  the 
given  initial  flight  condition.  First,  ground  collision  is  near 
immediate  eliminating  all  hope  for  ejection  success.  Next,  the 
uncontrolled  pitch  and  yaw  dynamics  illustrate  the  highly  unstable 
oscillatory  modes  that  are  to  be  expected  in  the  absence  of  active 
attitude  control.  The  acceleration  radical  evaluation  illustrated  at 
the  bottom  of  the  figure  with  scalings  (+/-  1,+/-  2)  indicate  that  the 
lethal  limit  is  attained  Immediately  (radlcal>l)  with  no  hope  for 
recovery.  A  thorough  inspection  of  auxiliary  hybrid  outputs  not 
illustrated  in  Figure  7.  1  indicates  that  the  maximum  analog  values  for 
several  variables  are  exceeded  almost  immediately  leading  to  the  loss 
of  simulation  fidelity  in  generation  of  the  final  results.  It  should 
be  noted  that  the  acceleration  radical  evaluations  appear  less  than  0 
in  the  strip  chart  outputs  due  to  severe  violations  of  allowable 
variable  bounds  in  the  analog  circuitry.  An  all  digital  VAX 
simulation  of  this  scenario  for  the  uncontrolled  case  indicated  that 
the  maximum  value  of  the  acceleration  radical  exceeded  10. 

The  PE  8/32  simulated  performance  labelled  "S"  in  Figure  7.1 
illustrates  clearly  the  advantages  of  active  translation  and  attitude 
control  when  compared  to  the  uncontrolled  performance  labelled  "U". 
Ground  collision  is  clearly  avoided  given  the  monotonical ly  increasing 
altitude  history.  The  desired  0  degree  roll  and  yaw  angles  are 
tightly  regulated  Indicating  proper  attitude  control  and  resulting 
ejection  seat  rotational  stability.  The  pitch  euler  angle  history  is 
indicative  of  successful  pitch  angle  control  given  the  well  controlled 
history  from  the  initial  0  degree  value  to  the  desired  nominal 
terminal  value  of  10  degrees.  The  rapid  reaction  of  the  control 
system  to  the  environmental  translational  forces  and  rotational 
torques  are  evident  in  this  case  which  completely  avoids  all  the 
negative  aspects  of  the  uncontrolled  ejection  seat.  The  resulting 
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acceleration  radical  is  bounded  to  less  than  1  throughout  the  2  second 
ejection  duration  with  a  terminal  reduction  by  an  order  of  magnitude 
of  the  acceleration  radical  and  dynamic  pressure  by  mission  completion 
when  compared  with  the  starting  values. 

The  microprocessor  hardware  in  the  loop  operation  labelled  "H" 
results  in  a  close  match  with  the  simulated  PE  8/32  results  labelled 
"S",  Some  discrepancy  is  observable  in  the  acceleration  radical 
history  between  the  "H"  and  "S"  results.  A  close  examination  of  the 
euler  time  histories  at  a  finer  strip  chart  scale  than  that 


illustrated  in  Figure  7.3  Indicated  some  variation  particularly 
noticeable  in  the  yaw  channel.  The  primary  disturbance  responsible 
for  this  variation  is  the  ambient  electrical  noise  at  the  simulated 
sensor  ADC  interface.  The  PE  8/32  operates  with  an  interface  voltage 
of  +/-  100  V  range  for  all  process  variables  generated  by  the  analog 
computer  while  the  Motorola  68000  ADC  interface  is  bounded  by  a  +/-  1C 
V  range.  The  RMS  electrical  noise  level  at  the  electrical  interfaces 
are  comparable  for  both  systems  Implying  an  order  of  magnitude  poorer 
signal  to  noise  ratio  for  the  Motorola  68000  sensor  measurements  when 
compared  to  those  available  to  the  PE  8/32. 

7.3  Low  Altitude  Escape  Conditions. 


The  initial  ejection  seat  flight  conditions  forming  the 
MIL-S-9479B  test  cases  are  defined  in  Table  7. 1.  The  corresponding 
figure  identifiers  illustrating  the  results  of  each  case  are  also 
indicated  in  the  figure. 
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TABLE  7.  1 

Low  Altitude  Escape  Conditions  (1) 


MIL 

Attitude 

Speed 

Altitude 

Results 

Case 

Pitch 

Roll 

Illustratii 

(Deg) 

(Deg) 

(KEAS) 

(.  feet ) 

1 

0 

60 

120 

0  (2) 

Figure 

7.2 

2 

0 

180 

150 

200 

Figure 

7.3 

3 

0 

0 

150 

300  (3) 

Figure 

7.4 

A 

-60 

0 

200 

500 

Figure 

7.5 

5 

-30 

0 

450 

500 

Figure 

7.6 

6 

-60 

60 

200 

550 

Figure 

7.7 

7 

-A  5 

180 

250 

600 

F igure 

7.8 

(1) 

Conditions 

at  start  of  escape  sequence 

(2) 

Impact  occurs  at  instant  of 

seat/aircraft 

separation 

(3) 

10000  feet 

per  minute  sink 

rate. 

All  the  cases  indicated  above  are  likely  to  result  in  ground 
collisions  if  the  opportunity  for  immediate  corrective  control  action 
is  missed.  Particularly  pressing  are  MIL  Case  3  with  a  high  sink  rate 
and  MIL  Case  7  with  inverted  roll  and  steep  pitch  orientated  to  cause 
ground  impact  in  the  absence  of  direct  control.  All  the  above  cases 
are  essentially  at  low  dynamic  pressure  with  the  exception  of  MIL  Case 
5  with  moderate  dynamic  pressure. 

Inspection  of  the  results  of  the  following  pages  for  the  "U"  or 
the  "uncontrolled"  class  of  response  dynamics  demonstrates 
unacceptable  performance  for  a  variety  of  reasons.  In  all  cases,  the 
attitude  is  highly  unstable  leading  to  high  attitude  rates  with 
corresponding  effects  observable  as  oscillations  in  the  acceleration 
radical.  The  altitude  performance  is  unpredictable  with  cases  2,  5,  7, 
leading  to  obvious  ground  collisions.  Ground  collision  avoidance  in 


other  cases  are  effectively  accidental.  In  all  instances  the  mean 
acceleration  radical  is  an  increasing  function  of  time  indicative  of  a 
poorly  controlled  situation. 

The  "S"  or  simulated  hardware  (PE  8/32)  results  by  contrast  have 
well  controlled  attitude  profiles  and  in  all  instances  ground 
collisions  arc  clearly  avoided.  The  terminal  pitch  angle  constraint 
is  10  degrees  for  all  cases  while  the  terminal  roll  angle  is  specified 
to  be  0.  In  all  cases  the  initial  yaw  angle  is  0  degrees  and  the 
intent  is  to  cightiy  maintain  0  degrees  yaw.  Inspection  of  the  MIL 
Case  7  results  illustrates  the  simultaneous  coordination  of  the  roll 
and  pitch  maneuvers  through  large  angle  rotations  with  minimal 
oscillations  demonstrating  the  effectiveness  of  decoupling  the 
attitude  control.  The  success  of  large  angle  euler  steering  is  also 
demonstrated  in  MIL  Case  2.  It  is  evident  that  the  terminal  euler 
constraints  are  successfully  met  satisfying  that  performance 
objective.  The  acceleration  radical  is  maintained  below  1  in  all 
cases  and  in  general  is  minimal  at  the  termination  of  each  mission 
scenario.  The  acceleration  radical  growth  during  the  mid  course  of 
each  mission  is  induced  by  the  actuation  system  control  maneuvers 
necessary  to  avoid  ground  collisions  plus  the  trajectory  conditioning 
maneuvers  necessary  to  meet  the  terminal  attitude  constraints. 

The  "H"  or  68000  hardware  based  controller  results  in  all  MIL 
Cases  (1-7)  are  essentially  identical  to  the  simulated  hardware  test 
cases  (PE  8/32).  These  tests  are  seen  as  conclusive  with  respect  to 
demonstrating  the  proper  operation  of  the  68000  based  controller 
algorithm  demonstrating  feasibility  of  hosting  the  control  concept  in 
currently  available  off  the  shelf  hardware.  The  results  of  the  next 
section  only  illustrate  the  simulated  hardware  performance  with  the  PE 
8/32  computer  since  all  testing  has  indicated  a  close  correspondence 
between  the  simulated  hardware  and  the  actual  microprocessor  hardware. 


The  remaining  results  exercise  the  control  strategy  at  various 
altitude,  attitude  and  initial  speeds.  Table  7.2  defines  the  initial 
ejection  seat  flight  conditions  for  these  test  cases.  The 
corresponding  figure  identifiers  illustrating  the  results  of  each  case 
are  also  indicated  in  the  figure.  Inspection  of  Table  7.2  indicates 
that  the  majority  of  flight  conditions  be  evaluated  at  a  number  of 
initial  orientations  in  order  to  test  controller  sensitivity  to 
initial  attitude  conditions.  The  dynamic  pressure  in  these  cases 
range  from  low  for  cases  (4-10),  moderate  for  cases  (11-15)  and  high 
for  cases  (16-20).  The  high  altitude  conditions  at  600  KEAS  for  cases 
(  16-20)  imply  actual  ground  speeds  of  approximately  2300  fps.  The 
high  dynamics  with  respect  to  absolute  inertial  velocity  in  cases 
(16-20)  leads  to  significant  coupling  of  the  attitude  and 
translational  controls  resulting  in  challenging  cases  to  deal  with 
when  a  single  set  of  parameters  for  all  possible  flight  conditions  to 
be  encountered  are  used  in  the  controller. 

The  results  from  the  uncontrolled  ejection  seat  in  this  series 
have  negative  conclusions  concerning  stability  and  containment  of  the 
acceleration  radical  similar  to  the  conclusions  of  the  previous  two 
sections.  Further  discussion  of  the  open  loop,  uncontrolled 
performance  is  not  warranted. 

The  "S"  or  simulated  hardware  performance  results  are  emphasized 
for  this  series  of  tests.  The  strip  chart  recordings  are  to  be  found 
on  the  ensuing  pages. 


TABLE  7.2 

Variable  Initial  Attitude  Conditions 


Attitude 

ch 

g) 


> 

) 


Speed 

Altitude 

Results 

Roll 

Illustration 

(Deg) 

(KEAS) 

( Feet ) 

0 

0 

0 

Figure  7.9 

0 

600 

0 

Figure  7.  10 

0 

700 

0 

Figure  7.11 

0 

150 

5000 

Figure  7 .  12 

120 

150 

5000 

Figure  7. 13 

180 

150 

5000 

Figure  7 .  14 

120 

150 

5000 

Figure  7.  15 

180 

150 

5000 

Figure  7.  16 

120 

150 

5000 

Figure  7. 17 

180 

150 

5000 

Figure  7.  18 

0 

400 

5000 

Figure  7.  19 

120 

400 

5000 

Figure  7.20 

180 

400 

5000 

Figure  7.21 

120 

400 

5000 

Figure  7.22 

180 

400 

5000 

Figure  7.23 

0 

600 

45000 

Figure  7.24 

120 

600 

45000 

Figure  7.25 

180 

600 

45000 

Figure  1 ,2b 

120 

600 

45000 

Figure  7.27 

180 

600 

4  5000 

Figure  7. 28 

Case  I  is  the  standstill  ejection  scenario  and  inspection  of  the 
results  illustrates  the  low  dynamic  pressure,  the  minimal  acceleration 
radical  plus  the  well  conditioned  attitude  time  histories.  The 
altitude  progresses  from  the  initial  0  value  to  a  terminal  value  of 
approximately  50  feet.  Cases  2-3  are  low  altitude,  high  dynamic 
pressure  conditions  with  conclusions  quite  similar  in  content  to  those 
of  paragraph  7.2.1. 

Cases  4-10  all  represent  low  dynamic  pressure  scenarios  with 
large  angle  variation  in  the  initial  attitude.  All  scenarios  are 
exercised  at  5000  feet  so  that  the  threat  of  ground  collision  is  moot 
in  these  instances.  The  terminal  pitch  angle  is  selected  to  be  10 
degrees  for  all  scenarios.  The  ejection  flight  conditions  do  not 
warrant  any  changes  in  roll  angles  to  avoid  ground  collisions  and  in 
these  situations  the  control  algorithm  will  regulate  the  roll  angle  to 
the  initial  values.  Examination  of  the  results  reveals  well 
conditioned  pitch  control  with  essentially  monotonic  transition  from 
the  initial  to  the  final  specified  value.  The  variation  in  the 
acceleration  radical  is  indicative  of  the  varying  control  energy 
requirements  required  to  induce  pitch  maneuvers  and  control  attack 
angle  in  attaining  the  terminal  constraint  values.  In  all  instances 
the  acceleration  radical  is  constained  within  the  minimal  injury 
region. 

Cases  11-15  exercise  a  subset  of  the  cases  4-10  at  the  indicated 
higher  speed  of  400  K.EAS  in  Table  7.2.  The  attitude  response  is 
similar  in  form  to  the  results  of  the  previous  paragraph  with  no 
instabilities  apparent  in  the  strip  charts  results.  Once  again  the 
roll  angles  are  clamped  to  the  initial  values  given  the  threat  of 
ground  collision  is  ruled  out  by  the  substantial  altitude.  The  higher 
control  energy  required  to  steer  the  angle  of  attack  to  the  terminal 
conditions  is  under  the  influence  of  the  higher  initial  speed 
condition. 
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Figure  7.12  VARIABLE  INITIAL  ATTITUDE  CONDITIONS 


U-UNCONTROLLED,  S=SinULATION  (PE  8/32).  H=HAPvDWARE  (MC  68000) 
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Figure  7,14  VARIABLE  INITIAL  ATTITUDE  CONDITIONS 
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Figure  7.15  VARIABLE  INITIAL  ATTITUDE  CONDITIONS 
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U-UNCONTROLLED,  S=SinULATION  (PE  8/32).  H=HAPvI!WARE  (MC  68000) 


The  final  series  of  test  cases  under  examination  are  cases  16-20 
which  are  high  altitude  and  high  dynamic  pressure  ejection  conditions. 
The  600  KEAS  ed  specification  translates  into  high  ground  speed 


upwards  of  2200  fps.  Hence,  rotation  of  the  velocity  vector  in  this 
area  of  the  flight  envelope  requires  very  high  energy  requirements  as 
evidenced  by  the  growth  in  the  acceleration  radical  for  cases  19  -  20. 
Case  results  results  in  an  excessive  radical  with  a  maximum  value  of 
approximately  1.2.  The  case  20  results  are  marginal  at  best.  The 
maximum  value  of  the  acceleration  radical  may  be  contained  in  these 
cases  by  a  more  judicious  choice  of  the  terminal  attitude  conditions 
imposed  on  the  problem.  A  very  careful  assessment  of  the  terminal 
angle  of  attack  condition  for  successful  ejection  is  warranted  to 
insure  that  the  control  energy  requirements  for  safe  ejection  are 
minimized.  The  extensive  testing  has  demonstrated  the  need  to  offer 
quick  rotational  stabilization  of  the  seat  frame  while  minimizing  the 
attitude  steering  to  insure  that  the  control  energy  does  not  itself 
become  the  major  source  of  pilot  injury. 

7. 5  Basic  Conclusions  on  the  Scenarios  Investigated. 

Examination  of  the  results  of  the  last  few  sections  indicates 
without  doubt  the  specific  benefits  in  the  use  of  advanced  control  for 
ejection  seats  in  improving  the  probability  of  survival  over  the  use 
of  simplified  booster  schemes.  The  method  of  control  summarized  in 
paragraph  6.3.1  yielded  superior  attitude  control  even  in  the 
instances  of  simultaneous  large  angle  corrective  steering  of  both 
pitch  and  roll  (see  MIL  Case  7  above).  The  high  altitude  response  of 
Cases  16-20  emphasized  the  need  to  select  the  terminal  attitude 
condition  carefully  in  order  to  minimize  the  control  energy  required 
in  attitude  steering.  The  low  altitude  performance  in  the  presence  of 
high  dynamic  pressure  (Cases  2-3)  were  wholly  acceptable  in  meeting 
the  control  objective.  From  a  control  standpoint  the  MIL  Cases  1-7 
are  the  most  challenging,  combining  the  effects  of  adverse  attitude 
with  imminence  of  ground  collision.  The  ability  of  the  control  scheme 
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in  dealing  with  those  eases  is  significant  evidence  of  the  viability 
of  the  method  in  handling  "difficult"  ejection  conditions. 


7,6  Microprocessor  Memory  and  Throughput. 

In  the  initial  program  development  plan,  the  PASCAL  language  was 
proposed  in  order  to  support  the  higher  order  language  (HOL) 
requirement  in  DOD  program  development.  An  HOL  environment  offers  the 
advantage  of  module  clarity,  ease  of  implementation  and  maintainance 
while  sacrificing  speed  of  execution.  The  HOL  strategy  leads  to  much 
lower  execution  rates  than  acceptable  so  that  assembly  code 
development  of  the  control  algorithm  is  required. 

Extensive  evaluation  of  the  execution  rate  requirements  for  the 
proposed  control  system  has  indicated  that  the  control  algorithm  need 
execute  at  a  minimum  of  a  40-50  Hz  rate.  Execution  rates  higher  than 
100  Hz  do  little  to  alter  system  performance.  The  results  of  the 
previous  section  were  all  generated  at  the  50  Hz  rate  demonstrating 
feasibility  of  operation  at  that  rate.  In  general  it  is  advisable  to 
dedicate  no  more  than  1/2  of  the  available  microprocessor  throughput 
to  control  algorithm  execution  leading  to  the  requirement  of  a  100  Hz 
algorithm  in  order  for  the  scheme  to  be  realizable  in  current 
hardware. 

Table  7.3  illustrates  the  memory  and  throughput  actually  measured 
using  the  TEKTRONIX  Microprocessor  Development  Lab  (MDL)  utilizing  the 
Motorola  68000  emulator  operating  with  an  8  MHz  clock  speed. 


TABLE  7.3 

Control  Algorithm  Throughput  and  Timing  Measurements  with 
Motorola  68000  Microprocessor.  (8  Mhz  System  Clock). 

Pascal  Code  Assembly  Code 

Memory  22  KBytes  20  KBytes 

Timing  0.28  s/cycle  (3.5  Hz)  0.007  s/cycles  (130  Hz)  ( 1) 

0.01  s/cycles  (100  Hz)  (2) 

( 1)  No  I/O  Included 

(2)  I/O  Included  (ADCs  &  DACs) 

These  timing  results  are  based  on  all  integer  arithmetic,  that  is, 
without  the  aid  of  the  faster  Motorola  68020  and  co-processor.  The 
original  PASCAL  code  results  indicate  a  nominal  3.5  Hz  execution  rate 
leading  to  the  conclusion  that  HOL  implementations  are  most  likely 
infeasible  on  most  available  state  of  the  art  hardware  components. 

An  assembly  code  version  of  the  control  software  was  developed  to 
evaluate  the  feasibility  of  control  algorithm  use  in  real  time 
applications.  The  measurements  indicated  in  Table  7,3  demonstrate  the 
capability  of  the  algorithm  to  execute  at  100  Hz  when  implemented  in 
assembly  language.  The  current  algorithm  is  then  feasible  given  50  Hz 
execution  rates  and  a  0.5  utilization  factor.  The  8  MHz  clock  speed 
is  intermediary  in  the  family  of  available  Motorola  68000  system 
'•locks.  Currently  the  Motorola  68000  may  be  utilized  in  conjunction 
with  any  of  the  following  clock  speeds:  4,  g,  16,  32  MHz.  Clearlv  the 
16  MHz  clock  would  not  push  the  system  clock  specification  to  the 
limit  and  would  make  feasible  200  Hz  algorithm  execution  speed  leading 
to  a  0.25  utilization  factor  of  the  processor.  The  16  MHz  clock  speed 
also  could  allow  for  alternate  tasks  to  be  executed  to  support  other 
functions  such  as  pyrotechnic  sequencing  and  drogue  chute  deployment. 
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Memory  requirements,  as  indicated  in  Table  7.3,  are  currently  in 
the  22  KByte  range.  These  memory  requirements  should  be  considered 
minimal  since  off  the  shelf  64  and  256  KByte  RAM  chips  are  readily 
available  at  low  cost.  Note  that  the  22  KByte  memory  requirement 
includes  the  total  local  and  global  data  storage  requirements  as  well 
as  the  controller  binary  code. 

In  summary  the  timing  and  memory  measurements  with  actual 
Motorola  68000  hardware  demonstrates  that  the  requirements  of  the 
control  approach  is  currently  feasible  at  relatively  low  expense  for 
the  computer  hardware.  The  low  accuracy  sensor  figures  of  Table  6.3 
of  the  previous  chapter  similarly  Illustrate  the  feasibility  of  use  of 
low  cost  inertial  components  in  supporting  the  hardware  requirements. 


7.7  Dynamic  Occupant  C.g.  Results 

In  compliance  with  para.  4.5.3  of  the  SOW,  SSI  has  Incorporated 
the  effects  of  pilot  motion  during  ejection  in  the  all  digital  VAX 
FORTRAN  simulation  in  order  to  reflect  c.g.  dynamics  of  the  combined 
seat/occupant  system.  Previous  work  aimed  at  demonstrating  the 
feasibility  of  attaining  a  desired  acceleration  tolerance  radical, 
attitude  and  altitude  given  the  assumption  of  a  fixed  c.g.  and  moment 
of  inertia  for  the  combined  seat/occupant  system.  In  reality  elastic 
properties  inherent  in  the  restraint  system  and  pilot  body  frame 
result  in  variations  in  the  c.g.  location  and  consequently  in  a  time 
varying  moment  of  inertia.  This  sub-task  is  intended  to  evaluate 
robustness  of  the  control  algorithm  in  maintaining  stability  of  the 
plant  and  containment  of  the  acceleration  radical  given  the 
disturbances  mentioned  above. 

Two  modifications  were  necessary  to  upgrade  the  existing 
simulation  in  order  to  account  for  the  dynamic  occupant  c.g.  The  first 
involved  the  implementation  of  the  seated  human  slump  model  forwarded 
by  R.J.  Dobbek  to  SSI.  That  model  represented  the  occupant  response 
to  forces  acting  on  the  ejection  seat  by  a  second  order  response  model 
in  each  of  the  three  orthogonal  directions  of  forward,  right  and  down. 
According  to  Dobbek,  the  original  model  was  obtained  from  the  SAFEST 
computer  program  documented  in  AFWAL  TR-82-3013  and  modified  to 
account  for  non-linear  characteristics  in  the  restraint  system  such  as 
the  "dead-zone",  "hysterises"  and  non  symmetrical  response  to  plus, 
minus  x  and  z  axis  loading.  Also  included  in  the  Dobbek  communication 
were  representative  values  for  the  parameters  in  the  second  order 
response  model  as  well  as  the  associated  non  linearities.  The 

variations  in  the  c.g.  dynamic  states  (positions  (6xcg,  Sycg.  iSzcg^* 

•  •  • 

velocities  (5xj,g,  ^Zcg^)  driven  by  the  second  order  response  model 

for  acceleration  (fix^g,  *5y2g,  Sz^-g)  were  appended  to  the  existing  6D0F 
flight  model  state  vector  consistent  with  the  Dobbek  model  and 
integrated  at  the  fundamental  integration  time  step  rate 
(approximately  0.001  s)  to  form  the  representation  of  pilot  c.g. 
motion. 
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The  second  modification  required  is  a  consequence  of  the  first, 
l.e.,  the  mass  shift  of  the  pilot  c.g.  results  in  a  time  varying 
moment  of  inertia.  Recall  that  the  center  of  gravity  of  n  masses  at 
time  t  is  expressed  by: 


_  mir^Ct) 

<'•» 

where 

™i  =  mass  of  body  i 

^i(C)  =  position  vector  for  center  of  mass  i  at  time  t  in 
some  reference  frame 

ri(t)  =  center  of  gravity  of  the  composite  of  the  system  in 
the  reference  direction 

The  second  moment  of  inertia  for  n  masses  about  an  axis  b  is  given  by: 

n 

2 


Ib(t)  =»  L  lai  +  mi(ri(t)  x  RCt))' 
i=l 

where 


(7.2) 


Ib(t)  =  moment  of  inertia  of  the  composite  of  n  bodies  about 
axis  b 

Ig,  ^  moment  of  inertia  for  body  i  with  respect  to  the  c.g. 
of  body  i  about  a  axis  a  parallel  to  b 

■^iC  t)  =  position  vector  of  body  i  in  some  reference  frame 

R( t)  =  unit  vector  defining  the  direction  of  axis  b  in  some 
reference  frame 


The  6D0F  dynamics  equations  previously  defined  in  are  valid 

when  the  moment  of  inertia  I  is  respect  to  the  system  c.g.  about  axes 
defined  along  the  orthogonal  directions  of  motion.  Hence,  R(t)  in 
equation  (7.2)  need  be  the  unit  vector  associated  with  the  center  of 
mass  as  defined  by  equation  (7.1)  for  Ib( t)  to  be  equivalent  to  I  in 
the  6  DOF  equations.  For  convenience  in  this  implementation,  ri(t), 
ri(t)  are  referenced  with  respect  to  the  SRP  frame.  The  only  time 
varying  term  in  the  expressions  (7.1,  7.2)  considered  here  is  the 
pilot  c.g.  location  in  SRP  units. 


The  data  for  evaluation  of  the  moment  of  inertia  of  the  combined 
seat/occupant  system  given  by  equation  (7.2)  is  taken  from  Air  Force 
supplied  data  illustrated  in  Table  7.4.  In  particular  the  evaluation 
here  assumes  the  95  percentile  pilot  data  from  row  8  while  the  data 
for  the  ejection  seat  is  given  by  the  composite  of  rows  1-7.  The 
initial  ejection  seat  flight  conditions  are  illustrated  in  Table  7.5 
so  that  the  dynamic  pressure  here  is  severe  implying  near  worst  case 
conditions. 

TABLE  7.5 

INITIAL  EJECTION  FLIGHT  CONDITIONS  FOR  DYNAMIC 
OCCUPANT  C.G.  EVALUATION 

Roll  0  degree 

Pitch  0  degree 

Yaw  0  degree 

Speed  687  KEAS 

Altitude  100  feet 

In  order  to  demonstrate  the  effects  of  the  occupant  c.g.  dynamics 
two  cases  for  preliminary  evaluation  were  constructed.  The  first 
assumes  the  pilot  c.g.  to  be  fixed  while  the  second  incorporates  the 
dynamic  effects  of  the  human  slump  model.  To  clarify,  the  distinction 
between  Cases  I  and  2  is  that  in  Case  1  the  system  c.g.  is  fixed  and 
the  control  law  has  an  exact  estimate  of  the  system  c.g.  and  moment  of 
inertia  while  in  Case  2  the  system  c.g.  and  moment  of  inertia  are  time 
varying  in  the  6  DOF  dynamics  equations  but  the  control  law  only  has  an 
estimate  of  the  initial  (an  approximate  time  average)  c.g.  and  moment 
of  inertia.  The  two  cases  then  directly  demonstrate  the  differential 
effects  of  the  additional  dynamics  due  to  pilot  motion  on  the  closed 
loop  control  system  performance. 

The  Case  1  results  are  shown  in  Figure  7.29  which  illustrate  the 
evolution  of  the  key  system  states  and  acceleration  radical. 

Similarly  the  results  from  comparison  Case  2  is  depicted  in  Figure 
7.30.  Inspection  of  the  two  cases  lead  to  the  preliminary  conclusion 
that  the  pilot  c.g.  dynamics  results  in  minor  variations  in  the 


control  system  performance  with  respect  to  control  of  the  plant  and 
acceleration  radical.  Perhaps  this  result  is  not  too  surprising  since 
the  maximum  excursion  of  the  system  c.g.  during  execution  in  Case  2  is 
bounded  by  2  inches  in  the  r.m.s  sense  with  a  maximum  deviation  of  1 
siug-ft2  in  the  moment  of  inertia.  Other  cases  with  comparable  level 
of  error  disturbances  of  the  same  order  demonstrated  a  similar  result. 


7.8  Nominal  Pilot  Results. 


This  section  examines  the  control  performance  in  an  all  digital 
simulation  environment  for  the  situation  when  the  controller  assumes 
the  parameters  associated  with  a  50  percentile  pilot  while  the  actual 
6  DOF equations  are  driven  by  parameters  associated  with  either  a  5  or 
95  percentile  pilot.  The  controller  only  has  approximate  knowledge  of 
the  seat/pilot  mass  and  moment  of  inertia  properties  which  is  the  case 
under  normal  conditions.  The  results  of  the  previous  paragraph 
focused  on  the  effects  of  the  CG  dynamics  due  to  pilot  motion  and  the 
conclusion  was  that  those  dynamics  were  Insufficient  to  alter  the 
basic  control  system  performance  measures.  The  cases  examined  here 
represent  large  constant  disturbances  imposed  on  the  control  system 
when  compared  to  pilot  motion  alone,  a  situation  of  potentially  far 
greater  consequences  with  respect  to  system  stability. 

All  previous  results  presented  so  far  were  based  on  a  limited 
definition  of  the  acceleration  radical.  The  results  presented  in  this 
section  incorporate  the  Dynamic  Response  Index  (DRI)  as  the 
acceleration  radical  so  that  a  more  realistic  evaluation  of  potential 
lethality  is  represented.  The  DRI  definition  and  associated 
parameters  are  based  On  the  requirements  established  as  part  of  the 
Crew  Escape  Systems  Technology  (CREST)  program.  In  that  program  the 
intent  is  to  produce  an  actual  prototype  of  the  next  generation 
ejection  seat  based  on  the  most  recent  advancements  in  applicable 
technologies. 

The  results  of  two  cases  are  presented  here.  Table  7.4  of  the 
last  section  defined  the  parameters  associated  with  the  ACES  II  seat 
on  a  component  basis.  Inspection  of  the  CG  shift  from  a  5  to  a  95 
percent  pilot  reveals  a  substantial  shift  for  the  z  CG  location.  The 
moment  of  inertia  components  for  the  respective  pilot  specifications 
indicate  a  substantial  difference  in  terms  of  the  x  and  y  principal 
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PILOT  RESULTS  FOR  NOMINAL  PILOT  DESIGN 


7.9  Robu>tn«i»;  Sensitivity  Analysis 


An  inportanc  aspect  of  any  control  system  is  Its  robustness. 

That  Is,  structural  parameter  changes,  sensor  failures  and  control 
component  defects  should  not  impede  the  controller  from  performing  Ics 
essential  tasVrs.  In  practice  it  is  desired  chat  Che  control  system 
remain  robust  against  a  given  set  of  failures  and  parameter  changes. 
Such  a  robustness  performance  is  usually  included  in  the  design 
characteristic  of  the  control  scheme. 

For  Vectored  Thrust  Digital  Control,  the  importance  of  its  robust 
behavior  stems  from  two  different,  although  related,  requirements. 
First,  the  Vectored  Thrust  Control  system  is  assumed  to  operate  under 
highly  uncertain  and  stressful  conditions  which  are  prevalent  at  the 
time  Chat  the  flight  crew  decides  to  eject  from  the  aircraft.  In  such 
conditions  parameters  change,  sensor  and  component  failures  are 
commonly  expected,  and  therefore  they  should  be  accounted  for  in  the 
design  phase.  Second,  one  operational  mode  of  the  controller  that 
has  been  described  in  previous  sections  relies  on  precomputed  and 
stored  trajectories  (reference  trajectories).  These  trajectories  and 
the  corresponding  control  sequences  are  computed  off  line  and  stored. 
Assuming  that  the  structure  of  the  seat  and  of  the  control  components 


remain  unchanged,  for  each  sec  of  Inlcial  condicions  (u,  v,  w,  p,  q, 
r,  X,  y,  8,9  ,  4i  ,i‘  )  one  crajeccory  and  its  associated  control 
sequences  have  to  be  stored.  It  is  clear  chat,  unless  some  data 
reduction  and  simplification  is  done,  Che  storage  and  retrieval 
problem  would  be  impractical  due  to  the  large  memory  requirement.  The 
data  reduction  and  simplification  has  to  rely  on  the  robustness  of  the 
control  system  with  respect  to  the  initial  condicions.  That  is, 
instead  of  computing  and  scoring  all  the  reference  trajectories,  a 
limited  number  of  them,  corresponding  to  typical  or  critical  initial 
points,  will  be  computed  and  scored,  together  with  the  corresponding 
control  sequences.  Then  for  any  inlcial  point,  a  reference  trajectory 
and  a  reference  control  sequence  will  be  approximated  (interpolated 
from  precomputed  trajectories  corresponding  to  the  neighboring 
points).  If  the  control  scheme  is  sufficiently  robust,  and  if  a 
sufficient  number  of  trajectories  have  been  stored,  it  is  expected 
chat  this  approximation,  coupled  with  an  on-line  regulation,  should 
result  in  a  reasonably  good  control  strategy.  In  the  following  we 
will  elaborate  on  various  aspects  of  this  issue. 

7.9.1  Derivation  of  Sensitivity  Relationships. 

In  this  section  we  derive  some  analytical  results  relating  the 
behavior  of  the  trajectory  to  the  change  in  initial  condicions  and  to 
the  change  of  the  control  vector. 

Let  us  define  the  following  notation: 

x(c)  :  crajeccory  vector  (12  dimenaiont) 

x(t)  -  (u,v,w,p,q,r,x,y,z,  9, 

f(.)  :  ^^12 

g(x)  :  ^gm 

U  :  m  dimenaional  control  vector 

Xq  :  initial  condition 

• 

*  ■  ^(x)  *  g(x)  U  :  nonlinear  dynamic  equation 
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x(t) 


0 


J 


f(x)dt  ♦ 


U(t)dt 


Robustness  with  respect  to  the  initial  vector  x^:  This  sensitivity 
matrix,  the  component  of  which  are  is  easily  computable  from 

the  trajectory  equation. 


Suppose  that  a  trajectory  x(t),  depending  on  the  control  vector 
U(t)  and  the  initial  condition  Xq  is  computed.  Assuming  that  Che 
control  vector  remains  the  same,  the  change  of  x(c)  relative  to  XqIs 
easy  to  compute.  In  fact: 

-  I+yi^Cxjdc  (x)U(t)dt  -  Si 

°  0  0 


vhcrs-|~  (x)  (x)  denote  ■eclces  of  derlveclve  of  the  elcaents  of 

f(x)  and  g(x)  with  respect  to  x.  Both  are  available  through  the 
linearization  of  the  dynamic  equation,  or  can  be  evaluated  on  line. 


For  example,  Che  above  formula  can  be  used  to  determine  the 
sensitivity  of  the  final  altitude  Z(T)  with  respect  to  different 
elements  of  Che  initial  conditions  x^.  This  altitude  is  independent 
from  Che  control  and  its  time  derivative  is; 


t  -  -sin  6  u  cos  6  sin  (^v  4-  cos  9  cos$  .w 


Then  it  follows 


^  Z(T)  f  , 

~  •  I  l-coiQ.u  -  ilne  ilni^.v  -  co«e  co«^.wldt 

0  J 


(co«9  co«$  V  -  co*0  tln^.wldt 


•ln9.dt 


az(T) 


cot  e  tin^  dt. 


cot  6  cot$  dc. 


Note  that  the  sensitivities  of  Z(T)  with  respect  to  the  initial 

vtlocltltt  dtptnd  on  cht  angles  0,  $  and  ip  rather  than  on  the 

■agnltudee  u^,  v  ,  w  . 

0  0  0 

In  general,  a  function  of  the  states  x  has  a  sensitivity  with  xq  which 
13  itself  a  function  of  the  elements  of  SI.  For  example,  RAD  is  a  function 
of  x  via  the  relationship 

RAD2  -  Ky  .  lfac(l)2  +  .  lfac<2)2  +  K„  •  lfac(3)2.  (7.4) 


ooo 


WWWW 
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In  (7.4  ),  Ifaci  Is  an  elsBanC  of  the  3-D  load  factor  vector,  Ifac.  This 
quancicy  is  defined  in  che  SAFEST  code  by  the  vector  equation 

Ifac  »  u)x(u)x  r^^g)  -  (u,v,w)  -  uxv  +  o)  x  tqq  (7.5) 

where  the  notation  in  (7.5)  is  standard.  The  in  (7.4)  are  constants  which 
define  nominal  load  factor  limits  for  pilot  survival  (we  have  used  a  con- 

•  e  • 

servatlve  value  for  Ky  in  place  of  the  DRI);  u,  v,  and  w,  in  (7.5)  are  ele¬ 
ments  of  f(x,u),  and  r^Q  is  a  constant  vector  in  this  analysis.  Equations 
(7.4)  and  (7.5)  are  then  used  in  the  RAD  sensitivity  matrix 


S3  -  —  (Ku 


RAD 


Ifac(I)  •  3lfac(l)/3x  +  K„  •  lfac(2)  • 


(7.6) 


31fac(2)/3x  +  Ky  •  Ifac(3)  •  31fac(3)/3x] 

In  deriving  (31fac/3x),  a  3-by-l2  matrix,  it  is  best  to  expand  each  of  the 

four  terms  in  (7.5),  and  then  take  explicit  partlals.  One  term  is  readily 

•  •  • 

available:  -3(u,  v,  w)/3x  is  the  negative  of  the  first  three  rows  of  the 
system  matrix  A.  The  other  terms  are: 
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where  9q/9x  is  the  fourth  row  of  the  system  matrix  A,  etc;  and, 


I Vdttf  /■*! 


^  (OJ  X  0)  X  rcc)  “ 


1  (qycG 

(pycG  ■  2qxcG) 

(PZCG  -  2rxcG) 

1 

03x3  1  ('P'CG  "  2pycG) 

(pxcG  +  rzcG> 

(qzCG  “  2r  ycG) 

,  03x6 

(rxcG  -  2pzcg) 

(rycG  -  2qzcG) 

(<iycG  p^^cg) 

1 

1 

where  (xcp,  y^c,  zqq)  are  components  of  rcg* 

It  Is  important  to  analyze  the  robustness  of  the  control  system 
with  respect  to  the  control  input  U(t).  In  practice  it  is  desirable 
to  use  the  same  pre-stored  control  input  U(t)  for  a  class  of  system 
conditions,  so  that  the  storage  and  retrieval  is  reduced.  On  the 
other  hand,  because  of  possible  defects  and  failures,  the  pre-stored 
control  vector  might  not  be  exactly  reallzeable.  If  the  system  is 
sufficiently  robust  a  crude  realization  of  the  control  input  U(t) 
should  result  in  an  acceptable  (although  degraded)  performance. 


Formally,  the  sensitivites  of  the  trajectory  x(c)  with  respect  to 
the  control  vector  U  can  be  calculated  from  the  nonlinear  first  order 
matrix  differential  equation: 


X 


U 


11 

3f 


Xy  U  +  g(X) 


Or  equivalently  from  the  metrix  Integral  equation: 

c  t  t 

*0*  y’g(x)dc  -  S2 

0  0  0 


where 


S2  -  X,, 


:  n  X  m  sensitivity  matrix 


The  above  CQueclons  can  be  solved  numerically*  In  most  cases  an 
approximate  value  of  the  sensitivity  factors  Is  sufficient,  and  can  be 
easily  obtained  either  by  numerical  methods  or  by  simulation. 

Our  experience  with  the  control  solutions  has  shown  that  the  control 
law  is  very  robust  with  respect  to  both  Initial  conditions  and  control 
input.  We  now  discuss  a  typical  example. 

Let  us  consider  a  typical  example  of  a  seat  trajectory  during  the 
course  of  the  ejection.  For  a  three  degree  of  freedom  model,  the 
optimal  solution  in  the  sense  of  maximizing  the  altitude  of  the  seat 
at  the  end  of  the  first  phase  of  the  ejection  is  shown  on  the  plots  of 
Figure  7.33  (a.b,c)  to  Figure  7.36.  The  state  vector  and  the 
corresponding  plot  symbols  are  defined  as  follows: 

Plot  Symbol 


u:  velocity  along  the  x  axis  XVEL 

w:  velocity  along  the  z  axis  ZVEL 

0  ;  pitch  angle  PITCH 

x:  position  along  x  axis  XFOS 

z:  position  along  z  axis  ZPOS 


The  control  vector  is  two  dimensional  and  is  defined  as 


U  ^  :  Thrust  level  (magnitude)  THRUST 

:  Nozzle  Angle 


NZLAMG 


The  iaecencaneous  DRl  Is  plocced  on  Figure  35a  ((denoted  by  DOBJ) 
and  Ice  Incegrel  (OBJFCN)  Is  plotted  on  Figure  7.36.  The  DRI  Is 
conscreined  Co  reaein  below  the  value  of  1,  alchough  in  practice  the 
violation  of  this  Halt  for  short  period  of  time  Is  tolerated. 

Using  the  above  solution  as  Che  base  case  trajectory  for  high 
speed  ejection  we  will  study  Its  variation  In  response  to  changes  of 
initial  conditions  and  of  Inputs. 

Test 

The  control  Inputs  remain  the  same.  The  Initial  velocities  u  and 
w  are  decreased  by  lOZ,  which  correpsonds  to  a  decrease  of  100 
f/s.  Figures  7.37  (a, b,c),  7.38  (a, b,c),  7.39  (a, b,c)  and  7.40 
correspond  to  this  test  case.  It  Is  clearly  seen  chat  the  pitch 
race  is  very  sensitive  to  these  initial  velocities.  The  two  most 
important  functions  to  monitor,  howerver,  are  the  final  altitude 
z(T)  and  the  DRI.  The  altitude  z(T)  decreases,  in  absolute 
value,  from  15  to  13  feet;  that  is  a  decrease  of  13Z  in  response 
CO  lOZ  In  Initial  velocities.  The  final  altitude  of  13  feet  Is 
acceptable.  The  DRI  function  change  Is  on  the  other  hand  very 
significant.  A  decrease  of  23  to  30Z  is  apparent  on  Che 
instantaneous  values  of  Che  DRI,  and  a  30Z  decrease  In  achieved 
on  Che  accumulated  (Integrated  DRI)  value  (Figure  7.38a  and  Figure 
7.40).  The  results  are  reasonable  and  agree  with  Intuition.  In 
fact  Che  thrust  level  being  maintained  at  its  previous  while  the 
initial  speed  Is  decreased,  it  is  natural  to  expect  the  seat  to 
decelerate  more  rapidly,  i.e.,  the  DRI  decrease.  The  same  Is  not 
time  If  the  velocities  were  to  Increase.  This  Is  illustrated  in 
Che  next  test  case. 

Test  2 

The  control  Inputs  remain  the  same.  The  Initial  velocities  are 
Increased  by  lOZ,  which  corresponds  to  a  100  f/s  increase. 


plots  of  Figures  7.41,  7.42,  7.43  and  7.44  illustrate  the 
trajectories  corresponding  to  these  test  case.  It  can  be  observed 
that  the  final  altitude  z(T)  Is  Increased  in  absolute  value,  by  14%. 
This  Is  natural  since  the  altitude  depends  directly  on  the  Initial 
velocity  in  z  direction.  On  the  other  hand  the  DRI  function  Increases 
by  30%  resulting  in  an  unacceptably  high  level  of  pressure  on  the 
pilot.  In  fact  since  the  thrust  level  has  not  been  augmented  to 
compensate  for  the  change  in  velocity,  the  acceleration  radial  remains 
relatively  high  on  the  pilot.  These  results  combined  with  the  ones  of 
the  previous  run  suggest  that  to  make  the  control  performance  robust 
against  increase  of  initial  velocity,  one  has  to  apply  more  thrust 
level  than  what  is  exactly  needed  by  the  ESOP  program. 

Test  ^ 

The  initial  conditions  remains  unchanged.  The  thrust  level  is 

smoothed. 

As  plots  of  Figure  7.33  display,  the  optimal  inputs  (optimal 
thrust  and  nozzle  angle)  computed  by  the  ESOP  programs  are  not  smooth. 
It  is  important  to  study  the  behaviour  of  the  solution  when  a  rather 
crude  smoothing  of  these  Inputs  is  performed. 

Figures  7.45  through  7.48  correspond  to  a  test  in  which 
the  thrust  level  la  smoothed  (see  Figure  7.46b).  The  essential 
features  of  the  trajectory  are  not  degraded.  In  particular  the 
altitude  z(T)  is  decreasd  by  I  ft.  which  corresponds  to  7%  decrease 
in  absolute  value.  The  performance  of  the  DRI  is  slightly  superior  to 
the  base  case.  It  is  interesting  to  notice  that  the  discontlnuties  of 
the  DRI  are  closely  correlated  with  those  of  the  thrust  level  afnd  the 


Test 

The  initial  conditions  remain  unchanged.  Both  the  thrust  level 
and  the  nozzle  angle  are  smoothed. 


Figures  7.49,  7.50,  7.51  and  7.52  correspond  to  the  case  where 
both  the  thrust  level  and  the  nozzle  angle  are  smoothed  (Figure 
7.49b,c).  As  expected  the  DRI  performance  remains  almost  similar  to 
the  previous  case,  with  somehow  less  dlscontinutltles.  The  final 
altitude  z(T)  Is  not  significantly  altered  by  the  smoothing  of  the 
nozzle  angle.  The  overall  trajectory  remains  accpetable  and  the  seat 
can  be  safely  ejected  using  these  two  smoothed  controls  Instead  of  the 
ones  generated  by  the  ESOP  program. 

7 .9 .2 .  Some  Implementation  Details.  It  Is  noted  chat  Che  expressions 
for  SI,  S2  and  S3  are  matrix  Integral  equations.  In  a  linear  system, 
closed  form  solutions  exist  or,  equivalently,  frequency  domain  results 
are  obtainable  via  Che  Laplace  transform.  Previously,  we  have  been 
integrating  a  system  of  12  equations.  The  introduction  of  the 
equations  for  Che  S  greatly  expands  the  system  order  to 

12  +  12^  +  (12  X  8)  -  252 

(x)  (SI)  (S2) 

(S3  is  derived  algebraically).  In  order  to  maintain  a  basic  level  of 
efficiency,  the  SI  system  was  made  more  compact.  Several  elements  of 
SI  remain  Identically  0  or  1,  thereby  requiring  no  integration.  These 
can  be  eliminated  from  the  system.  Doing  this  reduces  the  overall 
dimension  from  252  to  175,  a  30X  savings. 

The  rest  of  the  mechanization  Is  straightforward;  the  new 
Integration  equations  depend  strongly  on  the  linear  system  matrices,  A 
and  B.  These  are  available  from  existing  SSI  software  (as  well  as 
from  EASIEST),  at  any  point  in  time  along  the  trajectory.  The 
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approprlace  iniclal  conditions  are,  for  SI,  the  12-by-12  identity 
matrix  I,  and  S2(0)  »  B  At;  S3(0)  is  derived  algebraically  once 

ocher  quantities  are  developed. 

The  expanded  sensitivity  system  runs  at  a  noticeably  slower  rate 
than  Che  basic  system,  but  can  still  be  run  effectively  interactively. 

The  same  IMSL  OGEAR  algorithm  used  before  is  still  good,  with  much  the 
same  parameter  values. 

7 .9 .3 .  Results .  The  high  Q  (dynamic  preasure)-low  altitude  case  and  MIL-S 
Case  I  (no  dive,  60  deg.  bank,  120  KEAS  at  ground  level)  were  used  as 
references.  The  state  trajectory  is  in  agreement  with  the  original  result 
to  within  U.  The  variation  is  felt  to  be  due  to  the  automatic  changes  in 
integration  step  size  generated  by  OGEAR.  Tables  7.7  to  7.10  show  typical 
resul ts . 

Table  7.7  shows  SI  for  the  high  Q  case  at  t  -  0.96  second.  Each  ele¬ 
ment  of  SI,  3x|/3xqj  ,  Indicates  the  sensitivity  of  (row  i  of  SI)  to  a 
unit  change  in  (column  j  of  SI),  at  time  t  (0.96  sec.).  Only  67  of  the 
144  SI  elements  are  actually  integrated,  because  the  others  remain 
constant,  as  noted  above;  however,  all  144  SI  elements  are  displayed  for 
convenience.  Note  also  that  columns  defined  by  angular  quantities  (e.g.  , 
q,  pitch  race)  typically  have  larger  values  than  the  other  columns.  This 
is  due  to  the  radian  measure  used  for  angles.  For  example,  the  change  in  y 
(row  10  of  SI)  due  to  a  one  degree  change  in  ip,  heading  angle  (column  12) 
is  about  15  feet  at  0.96  seconds. 

Table  7.8  shows  MIL-S  Case  1  at  0.25  second.  Wliile  this  is  a  dif¬ 
ferent  case,  it  is  generally  true  chat  the  elements  of  SI  increase  with 
time.  Table  7.9  shows  the  same  SI,  with  S3  now  added  to  the  printout. 
Analyzing  one  of  the  elements  of  S3,  we  note  tat  the  sensitivity  of  RAD  to 
a  unit  change  in  Tq  (yaw  rate)  is  -0.0943.  In  general,  this  is  reasonable, 
because  Che  slideslip  angle  is  positive,  and  thus  a  positive  yaw  rate  would 
point  the  seat's  x  axis  more  into  the  wind  direction,  thus  reducing  RAD. 


I'ablf  7.7.  g,  Ixjw  Altitude 


lOOiH 

Olltl 


Table  7.10  shows  the  MIL-S  Case  I  at  t  -  0.09  second,  with  S2  added  to 
the  printout.  The  columns  of  S2  correspond  to  the  8  control  vector  ele¬ 
ments.  Columns  3,  6  and  8  of  S2  are  much  smaller  than  the  others,  due 
again  to  the  fact  that  the  latter  are  nozzle  deflection  commands  (radians), 
and  that  columns  3,  6  and  8  are  thrust  magnitudes  (pounds). 

7.9.4.  Conclusions ,  The  above  analysis  and  example  reveals  that  the 
robustness  of  the  control  scheme  should  be  studied  and  taken  into 
account  in  the  design  level.  More  precisely,  the  robustness  of  the 
control  structure  with  respect  to  initial  conditions  can  be  used  both 
for  data  reduction  and  storage,  and  also  for  deducing  on  line 
approximative  control  adjustment.  These  adjustoMnt 
procedures  are  of  simple  regression  type  and  can  be  implemented  on 
line. 

The  robustness  of  the  control  structure  with  respect  to  the 
control  inputs  should  be  taken  into  consideration  in  conjunction  with 
Che  technological  feasibility  and  the  cost  of  generating  the  desired 
wontroi  sequences.  For  instance* a  rapidly  oscillating  thrust  level 
may  not  be  easily  realizable,  and  one  would  like  to  replace  it  with  a 
less  demanding  thrust  level  variations  without  vitally  imparting  the 
overall  control  performances.  The  robustness  of  the  control  strategy 
can  be  taken  advantage  of  in  investigating  and  defining  such  a  control 


3 


oooooooc>e  rvioo 
f~  C'  try  r>  r  r)r»r'^oi->r» 
oooooooo<o^oo 

•  ••••••••••• 

00000©©Ok'^.0^ 


Of  OOOOOOOOOO 
oooooooooooo 
o<^for»ooor»coo 
OOOOOOOOOOOO 

•  ••••«••••«• 

oooooooooo^o 


OOOOOOOOOOOO 
c*  O  c  O  C»  r»  C»  o  C>  w'  o 
OOOOOOOOOOOO 
r't*'r»r*oo"ior  oc'o 

0''>000o000*-0  0 


O'lOCTOOOOO  "STCi 
0*-'0©OOiJOO  30Cj 

OOOOCJCjOOOOUO 

C‘Oi..oc:>oOO^>  c*o 


tO<y 

>  »  ^  't  c  '  J  ^  ^  f 

»•* 

S'  /X  O  O  «  I  /'  «~1 


<  1  ►»  r»  cj  .>  t  •  -jT  r  V'  »1  ^ 
Uxj'j 

<,->  4  V  •  X.'  .  .  •>  X#  a  xi  .  . 

««'  <JF''000^  J 


JOW^OOOO*^ 
-»«^C'-^OOOx'>  ■> 


^  f4  r  ^  .•»'.»x»  «' 

v'O*^ 


o'-ooc*«'>o©ooo6 

Wi  JUO0OU90OO00 
t>  '>**‘OO--3OOO000 

-y  '''-»OOr>©0000 


0;>irx^  0v*000000 
x»^*>^.0•4◦OUOOO 

r>''^o»^oooooo 
fv  ^^of^ooeoooo 

C>xi^OOuOOOOOOO 


w»><^w^OOOO 

.  ^  •  .->  r«  o  O  3  o 

j  •>  >^“30000 
'  >/  f  y  y  “y  'i  cy  Q 


.  >  »  fx,  rv 

o  ^  -  -/  i> 

f  •  ^  .  rx  o  o  r> 

*>•-  %'C/OUOOO 

.  wOOXi^W  Xa/^O 


^X  ^  f.  j  ^  ->  o  o  o 
/  ^  >  >  .•  >  o  O 


./X-  QOiJOOOO 

.  -^  r->  O  O  '•»  o  O  O 

>_Ow>OOC  JOO 
•  JOOuOOOOO 


»  ^  ■‘*9’  33^393 

»  ^  .X  n  •>  :»  -j  o  3 

.  s^*-9O.>9;j90 
.  f*  —  .»  3.>''r>'DO 


»-  j  O  J  J  J 

/  ■  I  1  5  1 

I  i  J  i  ^  J  -J 


■  jy  fs  <~y  ■>o->'a 
•  ->  O  i3^3  900;.»Q 

.  o  •<  «A  o  r>  o  r>  v*j  .3 

I  ««‘>^OO.IO90 


j 

>  j  •  •*  t  t  t  t  y 

►  III  j 


%% 


XVEL  IXIO'  1  ZVEL  (X10‘  I  PITRflTlXlO 


2.00 


ilii; 


SlJ  tv. 


TABLE  1.  Low  ALTITUDE  ESCAPE  CONDITIONS  d) 


PITCH 

ATTITUDE 

ROLL 

speed 

altitude 

(deg) 

(deg) 

(knots) 

(feet) 

0 

60 

120 

0  (2) 

0 

180 

150 

200 

0 

0 

150 

300  (3) 

-60 

0 

200 

500 

1 

CD 

0 

450 

500 

-60 

60 

200 

550 

-45 

180 

250 

600 

CONDITIONS  AT  START  OF  ESCAPE  SEQUENCE 

IMPACT  OCCURS  AT  INSTANT  OF  SEAT/aI RCRAFT  SEPARATION 

lOCGO  FEET  PER  MINUTE  SINK  RATE 
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Figure  7.68  MIL  7  (Table  1) 
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Figure  7.72  Case  4 
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Figure  7.74  Case  6  (Table  2) 
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Figure  7.79  Case  11  (Table  2) 
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Figure  7.^^5  Case  17  (Table  2) 
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CHAPTER  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

We  now  present  the  main  conclusions  from  the  design  effort  presented  In 
this  report: 

8. 1  Conclusions. 


The  Ejection  Seat  vectored  thrust  control  problem  Is  highly  nonlinear  and 
cannot  be  solved  effectively  or  efficiently  by  standard  linear  and  most  modern 
control  design  methods.  Other  multivariable  design  approaches  of  merit  are 
worthy  of  consideration  for  augmenting  the  approach  presented  here.  All 
reasonable  approaches  are  based  on  the  same  underlying  principles  of  robustness, 
disturbance  rejection  and  predictive  behavior.  Also,  the  nonlinear  equations  of 
motion  have  been  derived,  and  an  analytical  derivation  of  the  linearized 
equations  has  been  accomplished  and  partially  verified,  so  that  a  practical  and 
comprehensive  ejection  seat  control  design  and  analysis  package  now  exists. 

Also,  specification  of  three  hierarchical  levels  of  the  ejection  seat  control 
problem,  and  determination  of  the  methods  and  software  required  for  optimal 
control  design  at  all  levels  has  been  completed. 

The  EASIEST  program  has  been  used  extensively  for  computing  flight 
trajectories,  linearizing  models,  and  performing  simple  feedback  control 
strategies.  In  addition,  much  attention  was  given  to  quantification  of 
propulsion  system  energy  requirements  (Volume  I,  Chapter  3),  and  our  analysis 
presented  here  shows  that  high  Q  conditions  will  dictate  the  energy  needed  by 
the  propulsion  system,  although  some  low  speed,  low  altitude  adverse  attitude 
cases  also  require  considerable  energy.  At  least  8500  pounds  of  total  force 
capability  will  be  needed,  and  very  close  to  this  number  in  pound-seconds  of 
impulse.  This  number  could  be  revised  upward  If  thrust  vectoring  falls  much 
below  600  deg/sec  in  rate.  Also,  the  results  are  very  dependent  on  required 
reduction  In  dynamic  pressure  on  the  seat.  If  the  required  pressure  for  chute 


deployment  is  much  less  than  600  Q,  then  there  would  again  have  to  be  an  upward 
revision. 

Performance  results  in  Volume  II,  Chapter  7  show  that  the  safe  recovery  of 
the  pilot  under  extreme  ejection  conditions  is  feasible  using  vectored  thrust 
control  for  an  ACES  II  type  of  seat.  This  conclusion  is  based  on  examination  of 
at  least  thirty  extreme  escape  conditions  using  the  real  time  hybrid  simulator 
designed  for  this  project  (Volume  II,  Chapter  6).  The  thrust  levels  required 
are  achievable  and  the  control  solutions  are  robust.  Further  Improvements  are 
possible  by  using  a  flexible  propulsion  sytem  with  additional  thrust  rockets  and 
vernier  jets,  which,  most  likely,  would  be  used  to  achieve  the  required  degree 
of  reliability.  An  example  of  such  a  system  has  been  presented  in  Volume  I, 
Section  3. 2. 2. 4. 

Our  numerical  studies  to  date  have  pointed  out  the  following: 

High  Q  ejections  will  almost  certainly  be  aided  by  some  type  of  deployable 
fin,  to  exploit  aerodynamic  energy  for  attitude  stabilization.  Such  surfaces 
would  not  have  to  be  movable,  but  they  should  reduce  propulsion  system  required 
thrust  capacity. 

The  control  algorithm  has  been  successfully  implemented  and  tested  in  a 
real  time,  hybrid  environment,  hosted  on  a  wire-wrapped  controller  designed 
by  SSI  and  fabricated  by  Unidynamlcs.  The  simulations  confirm  the  validity 
of  the  design,  and  also  the  need  to  operate  the  controller  at  a  minimum 
sample  rate  of  50  Hz.  In  addition,  the  microprocessor  architecture  should 
require  no  more  than  about  half  of  available  microprocessor  throughput  to 
be  dedicated  to  controller  processing.  Memory  requirements  are  about  20 
Kbytes  for  assembly  code,  a  bit  higher  for  Pascal  code.  Pascal  code  will 
not  execute  the  controller  in  real  time  on  the  Motorola  68000  chip.  These 
results  confirm  that  the  controller  processing  requirements  can  be  met  with 
current  hardware. 


A  final  noce:  Che  control  concept  presented  in  this  report  has  become  a 
part  of  Che  major  technology  follow-on  project,  the  Crew  Escape  Technologies 
(CREST)  Program,  Working  as  a  subcontractor  to  Boeing,  SSI  has  overseen  Che 
application  of  Its  control  approach  to  an  actual  test  bed  seat.  This  project  is 
still  underway,  but  the  preliminary  results  are  that  the  control  logic  performs 
very  capably  on  Che  BMAC  baseline  design,  under  somewhat  more  difficult 
performance  criteria  than  chose  In  use  for  the  project  described  here. 

8.2  Recommendations. 


Even  though  significant  progress  has  been  made  in  demonstrating  feasibility 
of  the  vectored  thrust  control  concept,  further  testing  and  refinements  are 
necessary.  These  are  not  expected  to  be  technical  barriers  in  any  sense,  since 
the  most  difficult  issues  have  already  been  addressed. 

The  major  recommendation  is  that  Che  breadboard  simulation  system  developed 
under  this  project  at  Martin  Orlando  be  improved  and  expanded  upon.  At  this 
time,  it  Is  unique  In  terms  of  simulation  realism  for  analyzing  ejection  seat 
performance.  The  design  of  this  system  allows  for  efficient  progression  from 
just  the  wire-wrapped  controller  as  the  breadboard  hardware  element  (the  current 
design),  to  breadboard  systems  with  more  control  system  hardware  components 
added.  The  MMOA  High  Force  Gain  Valve  is  an  example  of  such  a  component  whose 
Integration  into  the  system  can  be  readily  analyzed.  The  simulation  capability 
described  in  Appendix  F  is  able  to  provide  timely  and  cost  effective  support  of 
advanced  design  projects  such  as  the  CREST  program. 
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Closed  loop  eigenvalues: 

-7.459  E3;  -7.459  E3;  -3.822  E+2;  -3.710  E2; 
-1,476  E2;  -2.238. E2;  -7.258  +  j  7.039; 
-1.393  El  +  j  1.396  El;  -1.;  -1.43 
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Table  8-7 
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Closed  loop  eigenvalues: 

-7.460  E3;  -7.460  E3;  -3.740  E2;  -3.564  E+2;  -2.240  E2: 

-1.330  E2;  -7.558  +  j  7.311;  -1.387  El  +  J  1.390  El;  -1;  -1.318 
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